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In  ejector  refrigeration  systems,  the  performance  of  the  ejector  is  critical  to  the  performance,  capability,  size 
and  cost  of  the  whole  system.  Construction  of  mathematical  models  has  been  used  as  an  effective  method  for 
analyzing  the  performance  of  the  ejector  as  well  as  the  whole  refrigeration  system.  These  models  can  also  be 
used  to  guide  system  operation,  interpret  experimental  results  and  assist  in  system  design  and  optimization. 
The  overall  objective  of  this  paper  is  to  provide  a  review  of  various  researches  of  the  mathematical  model  on 
the  hydrodynamic  and  thermodynamic  character  within  the  ejector.  The  paper  first  briefly  describes  ejector 
including  fundamental  principle,  flowing  and  mixing  mechanism  and  the  method  of  model  establishment. 
Then  various  models  consisting  of  ideal  assumptions,  governing  equations,  auxiliary  conditions,  solution 
methods  and  main  results  are  presented.  The  models  can  be  classified  into  two  main  categories:  (i)  steady 
thermodynamic  models  which  can  be  further  subdivided  into  single-phase  flow  model  and  two-phase  flow 
model  and  (ii)  dynamic  models  which  are  also  subdivided  according  to  the  flowing  phases  considered.  It  has 
been  shown  that  the  dynamic  models  have  higher  prediction  precision  and  give  more  information 
compared  with  the  steady  thermodynamic  models.  In  addition,  the  simplified  empirical  and  semi-empirical 
models  based  on  measured  data  are  briefly  discussed.  This  review  is  useful  for  understanding  the  evolution 
process  and  the  current  status  of  the  mathematical  models  on  ejector  and  highlighting  the  key  aspects  of 
model  improvement  such  as  the  mixing  mechanism,  the  capture  of  the  shock  wave,  etc. 

©  2008  Elsevier  Ltd.  All  rights  reserved. 


Contents 


1.  Introduction .  1761 

2.  Problem  description .  1762 

2.1.  The  fundamental  principle .  1762 

2.2.  Flow  and  mixing  mechanism  in  ejector .  1762 

2.3.  Mathematical  model  establishment .  1763 

2.3.1.  General  governing  equations .  1763 

2.3.2.  Auxiliary  conditions  and  mathematical  methods .  1763 

3.  Thermodynamic  model .  1764 

3.1.  Single-phase  flow  model .  1764 

3.1.1.  Constant-pressure  mixing  model .  1764 

3.1.2.  Constant-area  mixing  model .  1768 

3.2.  Two-phase  flow  model .  1769 

4.  Dynamic  model .  1773 

4.1.  Single-phase  flow  model .  1773 

4.2.  Two-phase  flow  model .  1775 

4.2.1.  Mixture  model .  1775 

4.2.2.  Eulerian  model .  1776 

4.2.3.  Interface .  1777 

5.  Empirical/semi-empirical  model .  1777 

6.  Conclusion .  1778 

Acknowledgements .  1779 

References .  1779 


*  Corresponding  author.  Tel.:  +8621  3420  6548;  fax:  +8621  3420  6548. 
E-mail  address:  rzwang@sjtu.edu.cn  (R.Z.  Wang). 


1364-0321/$  -  see  front  matter  ©  2008  Elsevier  Ltd.  All  rights  reserved, 
doi:  1 0.1 01 6/j.rser.2008.09.032 


S.  He  et  al / Renewable  and  Sustainable  Energy  Reviews  13  (2009)  1760-1780 


1761 


Nomenclature 

A 

cross-section  area  (m2) 

c 

local  sonic  velocity  (m/s) 

Cp 

constant-pressure  heat  capacity  (kj/kg  K) 

D 

diameter  (m) 

E 

total  energy,  Eq.  (121) 

f 

friction  factor,  Eqs.  (42),  (45),  (46),  (92) 

h 

specific  enthalpy  (kj/kg) 

k 

turbulent  kinetic  energy,  Eqs.  (4),  (124) 

l 

length  (m),  Eqs.  (45),  (46),  (92) 

M 

Mach  number 

m 

mass  flow  rate  (kg/s) 

P 

pressure  (MPa) 

R 

condensation  rate,  Eqs.  (108),  (112) 

Re 

Reynolds  number,  Eq.  (42) 

Kg 

universal  gas  constant  (kj/kg  K) 

s 

entropy  (kj/kg) 

T 

temperature  (K) 

u 

velocity  (m/s) 

U,  V,  W 

velocities  in  various  directions  (m/s) 

X 

location  along  the  ejector  axis(m),  Eqs.  (4), 

(94-96);  quality  of  the  fluid  (kg/kg) 

z 

axial  coordinate  in  cylindrical  geometry  (m) 

Greek  letters 

8 

void  fraction,  Eqs.  (106),  (109),  (110),  (115),  (117); 
rate  of  turbulent  energy  dissipation 

r 

heat  capacity  ratio 

0 

dissipation  energy  (kj/m3  s),  Eq.  (136) 

r 

production  rate 

T) 

isentropic  efficiency 

ll 

dynamic  viscosity  (Pa  s) 

P 

density  (kg/m3) 

X 

shear  stress 

V 

specific  volume  (m3/kg) 

Q 

area  contraction  ratio 

CO 

mass  flow  ratio 

momentum  loss  factor,  Eq.  (92) 

Subscripts 

c 

condenser 

d 

diffuser 

e 

exit;  evaporator 

f 

saturated  liquid,  Eqs.  (32)  and  (33) 

g 

generator;  saturated  vapor,  Eqs.  (32),  (33); 
gas  phase 

i 

inlet,  location  along  ejector  axis 

is 

isentropic  process,  Eqs.  (43),  (44) 

1 

liquid  phase 

lg 

phase  change 

m 

mixture;  mixing  section 

n 

primary  nozzle 

nt 

nozzle  throat 

P 

primary  stream 

s 

secondary  stream 

sat 

saturation 

ss 

conditions  downstream  of  shock  wave 

sub  sub-cooling 

t  total  temperature/pressure 

v  vapor 

y  mixing  section 

1,2,  3,  4  locations  along  the  ejector 


1.  Introduction 

The  ejector  refrigeration  system  has  great  market  potential 
because  it  offers  some  remarkable  advantages  which  include:  (i)  it 
can  alleviate  environment  problems  by  using  low  grade  thermal 
energy  sources  such  as  solar  energy,  geothermal  energy  and  waste 
heat  to  drive  the  system  instead  of  high  grade  electric  energy, 
hence  it  can  reduce  C02  emissions  resulting  from  the  combustion 
of  fossil  fuels,  (ii)  it  is  simple  and  with  no-moving  parts,  which 
makes  it  noise-free,  reliable,  long  lifetime,  low  initial  and  running 
cost  as  well  as  require  practically  no  maintenance  and  (iii)  natural 
substances  such  as  water  can  be  utilized  as  working  fluids,  which 
have  zero  ozone  depletion  potential.  The  main  flaws  of  the  system 
are  low  performance  efficiency  and  rigid  requirements  on  the 
operation  parameters  which  limits  their  widespread  application. 
The  most  direct  way  to  solve  this  problem  is  to  investigate  the 
ejector,  which  is  the  most  critical  component  where  highly 
irreversible  transformations  take  place.  In  particular,  due  to  its 
importance  on  the  whole  system,  the  ejector  can  be  taken  as  the 
core  of  an  optimization  analysis  aiming  to  define  its  operational 
conditions  and  the  geometric  configuration  that  make  the  COP  of 
the  system  maximum  [1].  Some  review  papers  have  been 
published  to  summarize  the  research  efforts  and  achievements 
focused  on  ejector  performance  optimization.  Chunnanond  and 
Aphornratana  [2]  summarized  the  ejectors  and  their  applications 
in  refrigeration.  They  concluded  that  the  understanding  of  the 
ejector  theory  had  not  been  completely  cleared.  In  order  to 
enhance  the  efficiencies  and  reduce  the  cost  of  ejector  cooling 
systems,  efforts  made  by  several  researchers  have  been  summar¬ 
ized  in  yet  another  review  paper  [3].  The  review  by  Sun  and  Eames 
[4]  outlined  the  developments  in  mathematical  modeling  and 
design  of  jet  ejectors.  However,  the  review  considered  only  the 
thermodynamic  model  based  on  two  basic  approaches,  the  mixing 
of  the  primary  fluid  and  the  entrained  fluid  either  at  constant- 
pressure  or  at  constant-area. 

With  the  development  of  computer  hardware  and  numerical 
methodology,  mathematical  models  are  being  used  for  better 
understanding  of  the  compression  process,  system  design  and 
performance  evaluation  considering  the  hydrodynamic  perfor¬ 
mance,  condensation  and  two-phase  flow.  The  advantages  of  this 
method  are  that  it  takes  less  time  and  cost  than  experimental 
method  for  predicting  the  performance  of  an  ejector.  The  second 
point  is  that  mathematical  models  can  produce  large  volumes  of 
results  at  virtually  no  added  expense  and  it  is  very  convenient  to 
perform  parametric  research  and  optimization  analysis.  The  third 
reason  is  that  by  experimental  means,  some  parameters  are 
difficult  to  be  obtained.  The  last  but  not  the  least  reason  is  the 
fundamental  physics  of  ejector  used  in  various  fields  are  similar,  so 
that  other  installations  may  be  investigated  too  by  utilizing  the 
model.  Consequently,  constructing  valid  mathematical  models  of 
the  ejector  has  become  the  key  subject  of  many  studies.  Lots  of 
mathematical  models  have  been  constructed  and  employed  to 
analyze,  develop  and  design  ejectors.  This  article  examines  the 
progress  made  in  the  area  of  mathematical  modeling  on  the 
ejector,  and  summarizes  comprehensively  the  numerous  signifi¬ 
cant  works  that  has  been  done  on  modeling  the  ejector. 
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Based  on  the  conservation  principles  of  energy,  mass  and 
momentum,  the  hydrodynamic  and  thermodynamic  performance 
of  ejector  can  be  formulated  with  complex  mathematical  detail  as 
governing  equations  and  auxiliary  relations.  This  is  the  main  and 
universal  method  for  constructing  mathematical  model.  There  are 
also  some  simplified  models  for  ease  of  computation  such  as  semi- 
empirical  model  [5].  Moreover,  it  is  an  effective  way  to  get  the 
empirical  performance  fitted  correlation  using  large  amounts  of 
experimental  results.  But  its  application  is  restricted  in  the  specific 
range  of  experimental  working  conditions.  This  paper  focuses  on 
the  review  of  the  mainstream  mathematical  models.  Besides,  some 
empirical/semi-empirical  models  are  briefly  introduced. 

The  first  section  of  the  paper  details  model  establishment  basis, 
which  can  be  either  thermodynamics  or  dynamics  arranged 
according  to  the  flowing  phases.  For  each  model  reviewed,  attempt 
has  been  made  to  gather  the  assumptions,  governing  equations, 
auxiliary  relations  and  the  solution  techniques  for  ensuring  the 
integrality  of  the  model.  The  results  and  validation  of  some  models 
are  listed  where  comparison  between  them  are  needed. 

2.  Problem  description 

2.2.  The  fundamental  principle 

An  ejector  is  also  known  as  jet,  injector  or  jet  pump  in  different 
literatures.  The  main  components  of  an  ejector  include  a  primary 
nozzle  (also  named  as  motive  nozzle  in  some  literatures),  the 
suction  chamber,  the  mixing  chamber  (including  a  convergent 
chamber  if  available  and  a  constant-area  throat  tube  which  is 
named  as  secondary  throat  in  some  literatures)  and  the  diffuser 
(Fig.  1). 

The  primary  nozzle  usually  is  a  convergent-divergent  nozzle. 
As  the  high  pressure  fluid,  known  as  “a  primary  fluid”  or  “motive 
fluid”,  expands  and  accelerates  through  the  primary  nozzle,  it 
flows  out  with  supersonic  speed  to  create  a  very  low  pressure 
region  at  the  nozzle  exit  plane  and  subsequently  in  the  mixing 
chamber.  Hence,  a  pressure  difference  between  the  streams  at  the 
nozzle  exit  plane  and  the  secondary  fluid  inlet  is  established  and  “a 
secondary  fluid”  or  “entrained  fluid”  is  drawn  by  the  entrainment 
effect. 

Based  on  Munday  and  Bagster’s  theory  [6],  it  is  assumed  that 
the  primary  fluid  flows  out  without  mixing  with  the  secondary 
fluid  immediately.  It  expands  and  forms  a  converging  duct  for  the 
secondary  fluid.  At  some  cross-section  along  this  duct,  the  speed  of 
the  secondary  fluid  increases  to  sonic  value  and  chokes,  were 
named  as  “hypothetical  throat”  (or  “Active  throat”  or  “aero¬ 
dynamic  throat”  in  some  literatures).  Then  the  mixing  process 
begins  after  the  secondary  flow  chokes  (double-choking,  critical¬ 
mode  operation).  It  is  assumed  that  the  pressure  of  the  two 
streams  is  uniform  at  the  mixing  section.  Also,  this  mixing  causes 


the  velocity  of  the  primary  flow  to  be  stepped  down  while  the 
secondary  flow  is  accelerated.  By  the  end  of  the  mixing  chamber, 
the  two  streams  are  completely  mixed  and  the  static  pressure  is 
assumed  to  remain  constant  until  it  reaches  the  constant-area  tube 
section.  Due  to  a  high  pressure  region  downstream  of  the  mixing 
chamber’s  throat,  the  flow  undergoes  a  succession  of  normal  and/ 
or  oblique  shock  waves,  also  called  a  shock  train,  accompanied  by  a 
corresponding  pressure  rise.  This  shock  causes  a  major  compres¬ 
sion  effect  and  a  sudden  drop  in  the  flow  speed  from  supersonic  to 
subsonic.  A  further  compression  of  the  flow  is  achieved  as  it  is 
brought  to  stagnation  through  a  subsonic  diffuser. 

It  can  be  observed  that  due  to  the  special  geometrical  structure 
of  the  ejector,  the  flow  in  the  ejector  is  very  complex,  which 
involves  supersonic  flow,  shock  interactions,  and  turbulent  mixing 
of  two  streams  and  two-phase  flow  in  some  cases.  Hence  the  main 
endeavor  of  modeling  the  ejector  is  to  reasonably  explain  the 
complex  flow  and  mixing  mechanism. 

2.2.  Flow  and  mixing  mechanism  in  ejector 

As  mentioned  in  Section  2.1,  the  flow  in  the  ejector  can  be 
divided  into  three  sections: 

(1)  Nozzle  section:  A  primary  fluid  is  accelerated  to  supersonic 
speed  by  the  convergent-divergent  primary  nozzle,  which 
forms  low  pressure  region  at  the  nozzle  exit  plane  and 
produces  the  entrainment  effect  to  entrain  the  secondary  fluid. 
Since  the  velocity  is  very  high,  the  process  can  be  considered  as 
adiabatic  and  the  isentropic  relations  are  applicable. 

(2)  Mixing  section:  Two  feasible  mathematical  methods,  the 
constant-area  mixing  model  and  constant-pressure  mixing 
model  proposed  by  Keenan  et  al.  [7]  are  used  to  describe 
the  mixing  mechanism.  Also,  a  normal  shock  wave  may 
occur  if  a  supersonic  fluid  mixture  enters  the  diffuser.  In  this 
case,  a  sudden  reaction  in  the  mixture  velocity  and  a  rise 
in  the  pressure  take  place.  In  addition,  the  fluid  mixture 
easily  undergoes  phase  change  and  a  condensation  shock 
may  occur. 

(3)  Diffuser  section:  The  mixture  of  primary  and  secondary 
flows  passes  through  the  diffuser,  and  converts  kinetic 
energy  into  pressure  energy.  At  the  diffuser  exit,  the  velocity 
is  reduced  to  zero  and  the  pressure  is  lifted  high  enough  to 
cause  discharge. 

Therefore,  the  mathematical  description  of  the  flow  inside  the 
ejector  is  very  complex.  Besides  the  conservation  equations  of 
mass,  energy  and  momentum,  some  gas  dynamic  equations,  state 
equations,  isentropic  relations  as  well  as  some  appropriate 
assumptions  need  to  be  used  to  assist  in  the  description  of  the 
flow  and  mixing  in  the  ejector. 


Fig.  1.  Principle  structure  of  ejector. 
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2.3.  Mathematical  model  establishment 

The  main  steps  used  to  establish  mathematical  models  include: 
(i)  assumption  proposal,  (ii)  deriving  the  governing  equations 
based  on  the  mass,  energy  and  momentum  balances,  (iii)  providing 
auxiliary  relation  to  close  the  governing  equations,  and  (iv) 
adopting  proper  mathematic  approach  to  solve  the  model.  For 
clarity,  the  system  is  considered  to  be  one-dimensional  (for 
thermodynamic  models).  Explanations  of  the  main  terms  in  the 
governing  equations  are  presented  below. 

Some  ideal  assumptions  are  needed  because  the  choking,  shock 
and  mixing  of  the  two  streams  occurring  in  the  ejector  are  too 
complicated  to  understand  completely.  On  the  other  hand,  some 
factors  that  do  not  influence  the  flow  significantly  can  be 
neglected.  Then  the  complexity  of  the  governing  equations  can 
be  reduced  and  the  solution  is  relatively  simple  and  time-saving. 
Some  basic  assumptions  made  before  establishing  mathematical 
models  are  as  follows: 

(1)  The  inner  wall  of  the  ejector  is  adiabatic. 

(2)  The  flow  inside  the  ejector  is  steady  and  isentropic. 

(3)  The  primary  fluid  and  secondary  fluid  are  supplied  to  the 
ejector  at  zero  velocity. 

(4)  The  velocity  at  the  ejector  outlet  is  neglected. 

(5)  The  two  fluids  begin  to  mix  with  a  uniform  pressure  at  the 
mixing  section. 

From  assumption  ( 1 ),  the  heat  transfer  between  the  ejector  and 
the  environment  can  be  neglected  and  the  isentropic  expressions 
can  be  used  based  on  assumption  (2).  Whereas  from  assumption 
(3),  the  inlet  fluid  is  at  a  stagnant  state,  therefore,  the  inlet  pressure 
and  temperature  of  the  fluid  are  also  equal  to  the  total  pressure  and 
temperature  as  well.  From  assumption  (4),  the  fluid  is  brought 
down  to  stagnation  when  it  discharges  out  of  the  ejector. 

2.3.1.  General  governing  equations 

Under  ideal  condition,  the  general  governing  equations  for  the 
nozzle,  mixing  and  diffuser  sections  are  as  follows: 

Conservation  of  mass  :  ^  P[U[A[  =  ^  peueAe  (1 ) 

Conservation  of  momentum  :  PXA[  +  ^  mxU[ 

=  PeAe  +  ^2  meUe  (2) 

Conservation  of  energy  :  y^  m^hj  +  uf  / 2) 

=  ^me(/i  e  +  ul/2)  (3) 


section  area  of  the  constant-area  throat  tube  are  also  used  for 
global  models.  In  addition,  the  inlet  conditions  at  the  boundary 
mesh  were  used  for  local  models. 

(3)  Turbulence  modeling:  As  the  flow  in  the  ejector  is  from 
subsonic  to  sonic  and  then  to  supersonic,  the  turbulence  has 
strong  effect  on  the  flowing  process.  Generally,  for  dynamic 
models,  the  Boussinesq  hypothesis  is  used  for  turbulence 
(Eq.  (4)).  This  means  that  they  are  based  on  an  eddy  viscosity 
assumption,  which  makes  the  Reynolds  stress  tensor  resulting 
from  equation  averaging  to  be  proportional  to  the  mean 
deformation  rate  tensor. 


■pu'^=^%+^ 


(/Jl c  -h  yut 


7ul. 


For  thermodynamic  models,  the  turbulence  cannot  be 
modeled  in  detail.  The  dissipation  term  is  interpreted  by 
frictional  and  mixing  losses  which  are  taken  into  account  by 
using  corresponding  coefficients  introduced  in  the  governing 
equations.  Generally,  these  coefficients  need  to  be  determined 
experimentally. 

(4)  Auxiliary  relations: 

Definition  of  Mach  number  :  M  =  u/c  (5) 


Definition  of  sonic  velocity  :  c 


dF 
dpj  js 

yRgT 


,  for  ideal  gas,  c 


(6) 


Based  on  isentropic  assumption,  the  following  gas  dynamic 
equations  are  frequently  used  in  the  thermodynamic  models: 


m  Pty/y[ 


A 


Mj 


1 


(Ki+l)/(Xi-l) 


M  v1  +ys1M? 


(7) 
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y+i 

2(y-l) 
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(9) 
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Entropy  equation  : 


Py 


constant 


(11) 


2.3.2.  Auxiliary  conditions  and  mathematical  methods 

In  order  to  close  the  equations,  some  auxiliary  conditions  have 
to  be  provided.  They  are  boundary  conditions,  initial  conditions  for 
iteration,  gas  dynamic  equations,  entropy  equation  and  the  state 
equation.  Also,  the  thermodynamic  and  transportation  properties 
of  the  fluid  should  be  given  as  well. 

(1)  Boundary  conditions:  Boundary  conditions  describe  the 
behavior  of  the  simulation  at  the  edges  of  the  simulation 
region.  The  boundary  conditions  commonly  encountered  in 
ejector  models,  are  usually  in  pressures  at  inlet  and  exit  of  the 
ejector.  Also  the  mass  flow  rates  or  the  velocities  of  the  primary 
and  secondary  fluids  have  been  used  as  boundary  conditions  in 
some  literatures. 

(2)  Initial  conditions:  The  initial  condition  used  in  ejector  models 
is  the  expansion  ratio.  The  entrainment  ratio  and/or  the  cross- 


State  equation  :  h  =  /i(P,  p)  (12) 

For  constant  pressure  mixing,  h  =  CPT  (13) 

(5)  Mathematical  solution:  Due  to  the  complexity  and  difficulty  in 
modeling  the  flow  in  the  ejector  completely,  most  of  the 
mathematical  models  are  one-dimensional  thermodynamic 
models  which  only  calculate  the  steady-state  explicit  equa¬ 
tions  and  obtains  state  and  operation  parameters  along  the 
ejector.  Detailed  information  such  as  shock  interactions, 
turbulent  mixing  of  two  streams  cannot  be  obtained  by  these 
models.  With  the  development  of  mathematical  method  and 
computer  science,  especially  with  the  birth  of  commercial  CFD 
software,  numerical  methods  have  been  applied  to  solve  the 
partial  implicit  differential  governing  equations  which  can  give 
more  information.  At  present,  the  most  widely  used  solution 
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method  for  ejector  model  is  finite  difference  method.  The 
governing  equations  are  solved  using  the  commercial  CFD 
software  package  FLUENT  with  the  governing  equations  being 
discretized  using  a  control  volume  technique. 


on  constant-pressure  mixing  theory  but  without  considering 
the  choking  of  the  secondary  flow.  For  simplification,  besides 
assumptions  (1-5)  mentioned  in  Section  2.3,  they  further 
assumed: 


3.  Thermodynamic  model 

Among  the  different  models,  the  model  based  on  thermo¬ 
dynamic  analysis  is  the  simplest,  which  solve  the  equations  in  one 
dimension.  In  this  model,  the  isentropic  relations  and  some  gas 
dynamic  equations  are  used  to  assist  in  the  description  of  the 
relationship  between  the  temperature,  pressure,  enthalpy  and 
Mach  number  and/or  velocity.  Fience,  most  of  them  are  explicit 
steady-state  models  and  work  at  design  conditions. 

Moreover,  in  order  to  simplify  modeling,  most  of  the  models 
assume  that  the  inlet  fluid  is  superheated  and  no  phase  change 
occurring  in  the  ejector.  These  models  are  developed  based  on 
single-phase  flow.  Flowever,  the  experimental  results  indicate  that 
the  fluid  in  the  ejector  easily  undergoes  condensation  shock.  Also, 
for  some  practical  application,  the  ejector  works  as  a  pump  by 
using  high  pressure  water  to  entrain  gas  or  high  pressure  vapor  to 
entrain  water.  In  these  cases,  the  flow  in  the  ejector  is  two-phase 
flow.  Therefore,  in  the  next  section,  the  models  have  been 
subdivided  into  single-phase  flow  models  and  two-phase  flow 
models. 


(1)  Friction  losses  were  introduced  by  applying  isentropic 
efficiencies  to  the  primary  nozzle,  diffuser  and  mixing 
chamber. 

(2)  The  two  streams  mixed  at  the  primary  nozzle  exit  plane. 

(3)  Mixing  of  the  two  streams  is  complete  before  a  normal  shock 
wave  can  occur  at  the  end  of  the  mixing  chamber. 

(4)  The  working  fluid  is  ideal  gas. 

(5)  The  flow  is  one-dimensional. 

(6)  The  constant-pressure  mixing  occurred  at  the  inlet  of  the 
constant-area  section  (Fig.  2). 

Based  on  these  assumptions,  the  model  was  then  written  as: 
Primary  nozzle: 

l 

Energy  equation  :  =  rj„(hpX  -  hpl)  (15) 


(16) 


3.2.  Single-phase  flow  model 

3.1.1.  Constant-pressure  mixing  model 

Since  the  birth  of  the  first  steam  jet  refrigeration  system,  its 
advantages  have  made  it  to  enjoy  popularity.  In  1942,  the 
theoretical  prediction  of  the  ejector  progressed  further,  when 
Keenan  and  Neumann  [8]  established  a  one-dimensional  con¬ 
tinuity,  momentum  and  energy  equations  to  predict  the  perfor¬ 
mance  of  the  ejector.  Flowever,  the  difficulty  in  offering  a 
reasonable  analytical  solution  for  the  momentum  conservation 
equation  (Eq.  (14))  during  mixing,  made  the  use  some  experi¬ 
mental  coefficient  necessary  in  model. 

=  j)  AdP  =  mp(l  +  co)um  -  mpu2 p  -  msu2s  (14) 

Therefore,  in  Keenan  et  al.’s  following  work  [7],  they  introduced 
two  feasible  theoretical  methods  to  solve  the  problem,  the 
constant-pressure  mixing  model  and  the  constant-area  mixing 
model.  The  former  considers  that  dP  =  0  during  mixing  while  the 
latter  considers  dA  =  0. 

Keenan  et  al.’s  work  lays  a  foundation  for  one-dimensional 
design  theory  of  the  ejector.  They  pointed  out  that  a  constant- 
pressure  mixing  ejector  gives  better  performance  than  a  constant- 
area  mixing  ejector.  Therefore,  most  of  the  mathematical  models 
that  followed  are  based  on  the  constant-pressure  mixing  ejector. 
However,  Keenan  et  al.’s  initial  model  could  not  give  detailed 
information  on  the  choking  phenomena  that  is  easy  occurring  in 
the  supersonic  ejector  when  the  back  pressure  is  low. 

Munday  and  Bagster  [6]  further  developed  the  constant- 
pressure  mixing  model  by  assuming  that  the  primary  fluid  flows 
out  without  mixing  with  the  secondary  fluid  immediately  after 
discharging  from  the  nozzle  exit.  They  assumed  that  both  fluids 
mix  somewhere  at  the  downstream  of  the  nozzle  exit  (Fig.  1, 
section  y-y)  in  the  suction  chamber  and  a  “hypothetical  throat”  for 
entrained  fluid  is  formed  here. 

However,  neither  Keenan’s  nor  Munday’s  model  has  taken 
consideration  of  irreversibility  due  to  friction.  Eames  et  al.  [9] 
modified  Keenan’s  model  to  include  irreversibility  associated 
with  the  primary  nozzle,  mixing  chamber  and  diffuser  based 


(17) 


where  hpl  is  the  flow  enthalpy  at  the  nozzle  exit  for  isentropic 
expansion. 

Mixing  section:  Based  on  assumption  (6),  Pi  =  Pm,  Aa  =  Am,  and 
assumption  (5),  the  momentum  equation  (Eq.  (2))  between  section 
1-1  and  section  m-m  can  be  simplified  into  Eq.  (18)  with  rjm  as  the 
efficiency  for  the  whole  mixing  chamber. 


Pm(rnPIipl  +  msUsi)  =  (mp  +  rns)um  (18) 

Using  Eq.  (6),  Eq.  (18)  then  can  be  written  in  terms  of  the  Mach 
number: 


_  Mpi  +  coM*^  V Ts/Tp 
m  _  (ru+  1)(1  +coTs/Tp) 


(19) 


where  co  is  the  entrainment  ratio:  &>  =  ms/mp,  and 

M*  =  y^(T±lM2)/^l  +  ^M2j,  it  is  the  actual  mixture  velocity 

divided  by  the  sonic  velocity  of  the  mixture  at  critical  conditions. 

Shock  wave:  The  relations  between  the  Mach  number 
and  pressure  lift  upstream  and  downstream  of  the  shock 
wave  are 


(20) 


P3/P  m 


1  +yM] 

1  +  yM2m 


(21) 


Diffuser  section:  The  pressure  lift  in  the  diffuser  described  in 
Mach  number  including  r/d  as  efficiency  is 


P4/P3 


(y  -  i)Pd 


Mo  +  1 


y 

y-i 


(22) 


Give  the  temperature,  pressure  and  mass  flow-rate  of  the 
primary  and  secondary  fluids  and  assuming  a  value  of  Pi/Ps>t, 
calculating  Eqs.  (16),  (17),  (19-22)  consequently,  then  you  can 
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Fig.  2.  Ejector  model  used  by  Eames  et  al. 


obtain  the  exhaust  pressure  P4,  repeat  the  calculations  with  a  new 
value  of  P1/PSit  until  the  maximum  P4  is  obtained.  The  primary 
nozzle,  diffuser,  and  mixing  chamber  efficiencies  of  0.85,  0.85  and 
0.95,  respectively  were  used  for  the  analysis.  The  authors  verified 
the  theoretical  results  with  experimental  data  and  found  that 
choking  of  the  secondary  flow  in  the  mixing  chamber  of  the  ejector 
plays  an  important  role  in  the  system  performance.  Maximum  COP 
can  be  obtained  when  the  ejector  is  operated  at  its  critical  flow 
condition. 

Based  on  the  theory  proposed  by  Munday  and  Bagster  [6] 
and  using  the  model  developed  by  Eames  et  al.  [9],  Aly  et  al. 
[10]  presented  two  models.  One  was  developed  by  applying 
steady-state  equations  of  energy,  momentum  and  continuity  at 
the  nozzle,  diffuser  and  mixing  section  to  determine  the 
pressure  and  velocity  at  each  section  from  which  the  system 
performance  can  be  obtained.  Based  on  the  assumptions  (1-5) 
mentioned  in  Section  2.3  and  taking  into  account  of  the  losses 
due  to  wall  friction  in  the  nozzle,  and  those  in  the  mixing  and 
diffuser  sections,  the  mathematical  description  of  the  model  is 
as  follows: 

Primary  nozzle: 

“2“  =  Pn(^p,t  —  frpl )  (15) 

Mixing  section:  In  the  mixing  section,  both  the  momentum  and 
energy  balance  equations  were  used.  Comparing  with  Eq.  (18),  the 
momentum  term  of  the  secondary  fluid  at  left  side  is  omitted 
(Eq.  (23))  based  on  assumption  (3)  in  Section  2.3.  The  energy 
balance  is  shown  in  Eq.  (24). 

=  (mp  +  ms)u2  (23) 


(mp  +  iris)  h2  + 


mshs  +  mp 


(24) 


Shock  wave:  The  Eq.  (25)  similar  to  Eq.  (20)  was  derived  to 
calculate  the  mixture  velocity  after  the  shock  wave 


l'3  =  y(^T  +  M22)/(^TM22  +  l)  xc3  (25) 

By  applying  the  energy  and  continuity  balance  before 
and  after  shock  wave,  the  enthalpy  and  density  after  the  shock 
wave  were 


Diffuser  section:  The  energy  equation  taking  account  the 
isentropic  efficiency  for  the  diffuser  section  is 


(28) 


where  /i4>is  is  the  flow  enthalpy  for  the  isentropic  process. 

Assuming  a  certain  value  for  the  entrainment  ratio  and 
considering  the  isentropic  expansion  process  in  the  primary 
nozzle,  the  pressure  Pi  and  enthalpy  h\  at  the  nozzle  exit  can  be 
obtained.  Then  Eqs.  (15),  (23-28)  can  be  used  to  obtain  the  diffuser 
exit  enthalpy  h4.  The  calculation  is  repeated  with  a  new 
entrainment  ratio  until  the  required  exit  pressure  is  obtained. 

The  second  model  used  the  assumptions  in  the  first  model  and 
in  addition,  the  flow  inside  the  ejector  was  considered  as  ideal  gas 
with  constant  specific  heat  ratio  both  for  the  superheated  region 
and  the  wet  region  separately.  Following  Eames  et  al.’s  model  of  a 
steam-vapor  ejector  [9],  they  calculate  the  pressure  ratio  and 
Mach  number  of  the  flow  in  the  ejector  to  study  the  characteristics 
of  the  system. 

For  validity,  the  results  of  the  two  models  were  compared  with 
the  empirical  correlation  of  Power  [11]  and  good  agreement  was 
found.  Then  both  models  were  used  to  analysis  the  effect  of  the 
efficiencies  such  as  nozzle  and  diffuser  efficiencies,  as  well  as  the 
design  parameters  on  the  performance  of  the  ejector.  They  found 
that  the  nozzle  and  diffuser  efficiencies  have  a  significant  influence 
on  the  system  performance. 

In  order  to  take  account  of  the  chocking  of  secondary  fluid  that 
neither  Eames  et  al’s  nor  Aly  et  al’s  models  considered,  Huang  et  al. 
[12]  presented  a  critical-mode  (double-chocking)  model  based  on 
Munday  and  Bagster’s  theory  [6]  by  assuming  that  the  primary 
fluid  mixed  at  some  cross  section  (Fig.  2,  section y-y)  downstream 
of  the  nozzle  exit.  Hence,  the  model  is  more  complex,  besides  all 
equations  derived  by  Eames  et  al.,  the  gas  dynamic  relations 
(Eqs.  (6-10))  and  state  equation  Eq.  (13)  were  used  in  the 
calculation.  In  addition,  the  area  relationship  (Eq.  (29))  and  the 
energy  conservation  in  mixing  section  (Eq.  (30))  were  used  as  well. 


Apy  +  A 


sy 


(29) 


mp 


+  iris 


^CpTsy  + 


=  (mp  +  ms) 


^CpTm  + 


(30) 


,  ui  ui 

h2+f-f 

(26) 

u2 

P2^3 

(27) 

Eleven  different  ejectors  were  tested  to  determine  the 
coefficients  used  in  the  model  and  to  verify  the  theoretical 
analysis.  Then  the  performance  analysis  was  carried  out  to 
determine  the  entrainment  ratio  and  the  required  cross  sectional 
area  of  the  constant-area  throat  tube.  In  this  model,  the  initial 
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Fig.  3.  Ejector  and  the  Molier  chart  of  an  ammonia  ejector. 


conditions  were  the  throat  and  exit  area  of  the  primary  nozzle,  the 
stagnant  temperature  and  pressure  of  the  primary  and  secondary 
fluid  and  the  critical  back  pressure. 

It  should  be  pointed  out  that  the  previously  mentioned  models 
are  based  on  ideal  gas  assumption  which  does  not  reflect  the  actual 
process  occurring  in  the  ejector.  To  solve  this  problem,  Rogdakis 
and  Alexis  [13]  improved  the  model  proposed  by  Munday  and 
Bagster  [6]  by  using  the  thermodynamic  and  transportation 
properties  of  real  gas.  The  model  was  based  on  the  Molier  chart 
of  an  ammonia  ejector  (Fig.  3(b)).  The  momentum  conservation 
during  mixing  was  employed  in  the  model.  Detailed  governing 
equations  are  as  follows: 


Primary  nozzle: 

The  energy  balance  : 

u2aJ2  =  hi  -/jal 

(31) 

The  state  equation  : 

Sl  =  S(Tj,Pj)  =  Sas 

-^asSgas  (1  —  *^as)Sfas 

(32) 

and 

has  —  ^ashgas  +  (1  —  ^as)hfas 

(33) 

The  isentropic  efficiency  is  rjlal  =  (hi  -  hai)/(hi  -  has) 

(34) 

Similar  equations  were  developed  for  the  secondary  fluid. 
Mixing  section: 

Momentum  balance 

:  Ua  i  +  (OUa2  —  (1  + 

(35) 

Similarly,  by  employing  the  energy  balance  and  state  equations 
Eqs.  (31 )  and  (32),  the  enthalpy  hc,hd  and  the  specific  volume  vc,  vd 
can  be  obtained. 


Diffuser  section: 

sd  =s(Td,Pd) 

(36) 

S3s  =  Sd  =  S(T3s,Pd) 

(37) 

h3s  =  h(T3s,  Pd) 

(38) 

nd3  =  (hd  -  h3s)/(hd  -  h3) 

(39) 

h3  =  h(T3,P3) 

(40) 

Given  the  pressure  for  the  state  (1-3)  (Fig.  3(b))  and  using 
iterative  calculation,  the  maximum  entrainment  ratio  oo  can  be 
found  with  constant  superheated  temperature  TV 

The  ejector  efficiency  of  compression  is: 

r\  =  (mg  +  me)(h3  -  h2)/(mg(hi  -  has)) 

=  (1  +  co)(h3  -  h2)/(hi  -  /ias)  (41 ) 

The  subscripts  1,  2,  3,  b,  c,  d,  al,  a2,  as  and  3s  are  the  cycle 
locations  as  illustrated  in  Fig.  3. 

The  model  can  be  used  for  one  component  or  two  components 
ejector.  They  used  the  model  to  analyse  the  performance  of  an 
ejector  cooling  system  providing  air  conditioning  and  it  was  found 
that  the  shock  phenomena  on  the  entrained  vapor  play  important 
role  in  ejector  performance  and  the  COP  of  the  system.  Also,  from 
their  work,  they  observed  that  the  highest  operating  efficiency  can 
be  achieved  if  the  generator  pressure  increases  with  decreasing 
condenser  temperature,  when  the  evaporator  temperature  is  fixed. 
In  addition,  they  showed  that  the  maximum  coefficient  of 
performance  is  a  linear  function  of  generator  temperature,  a 
quadratic  function  of  condenser  temperature  and  a  cubic  function 
of  evaporator  temperature. 
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Fig.  4.  Ejector  model  used  by  Selvaraju  and  Mani. 

The  friction  loss  is  closely  related  with  the  velocity,  and  the 
velocity  varies  in  a  big  range  along  the  ejector.  Therefore,  using  a 
constant  coefficient  to  represent  the  friction  loss  as  in  above- 
mentioned  models  is  a  simplified  treatment.  Also,  Aly  et  al.  [10] 
found  that  the  efficiencies  have  a  significant  influence  on  the 
system  performance.  Taking  this  into  account,  Selvaraju  and  Mani 
[14]  developed  a  model  based  on  Munday  and  Bagster’s  [6]  theory 
for  critical  performance  analysis  of  the  ejector  system  (Fig.  4)  by 
applying  an  expression  to  describe  the  friction  loss  in  the  constant- 
area  section: 

— !=  =  2.01og(J?em  /u  -  0.8  (42) 

V  J  m 


where  rjn,rjm  and  rjd  are  the  efficiencies  of  the  nozzle,  mixing 
chamber  and  the  diffuser.  They  can  be  obtained  from  the  following 
expressions: 


Pn 


frp,t  frpl 
^p,t  —  frpl,is 


(49) 


P  m  = 


u 


m 


U 


mp 


(50) 


Pd 


^md.s  —  7lmm 

fimd  —  fimm 


(51) 


This  model  was  used  to  predict  the  ejector  refrigeration  system 
with  mechanical  sub  cooling  which  uses  an  auxiliary  liquid-gas 
ejector  to  enhance  sub  cooling  for  the  refrigerant  from  condenser. 
It  was  found  that  the  performance  of  this  cycle  improved  more 
than  the  traditional  ejector  refrigeration  system  and  an  optimum 
mechanical  sub-cooling  degree  exists  for  different  refrigerant  and 
operating  conditions. 

In  order  to  obtain  more  detailed  information  on  the  ejector, 
Ouzzane  and  Aidoun  [16]  derived  a  local  mathematical  model 
and  computer  programs  for  ejector  studies  in  refrigeration 
cycles,  one  program  for  optimal  ejector  design  and  the  other  for 
simulation  with  more  in-built  flexibility.  The  model  is  based  on 
Munday  and  Bagster’s  theory  [6]  and  isentropic  flow  in  the 
nozzles  and  the  diffuser.  The  control  volume  along  the  ejector 
axis  is  shown  in  Fig.  5.  The  conservation  equations  for  the 
control  volume  are, 


Then  they  derived  the  minimum  area  per  mass  flow  rate  and  the 
velocity  for  both  primary  and  secondary  fluid  at  the  primary  nozzle 
outlet  and  the  aerodynamic  throat,  respectively. 


Vr 


A 

mn 


upe  —  y  2?7p(/lg  —  fipe). 

The  velocity  of  the  mixed  fluid  is  expressed  as, 


1  9  1  9 

Energy  conservation  :  hx  +  -uf  =  hx_ i  +  - uf_a 


Mass  conservation  :  p{U[A[  =  /oi_1  ui_1yAi_1 


(52) 


(53) 


At)* 

(43) 

Isentropic  condition  :  S\  = 

(54) 

hpe)is 

(44) 

Momentum  conservation  :  pi uf  +Pj  =  +Pi_i 

(55) 

The  iteration  variation  for  the  solution  procedure  is  Mach 
number,  such  that 


um 


(rripUpe  +  msume)  +  (P  se 


mm 


Pm  )Am 


(45)  Mj  =  Mj_!  +  AM 


(56) 


The  entrainment  ratio  and  the  characteristic  area  ratio  for 
ejector  are, 


Zipe  —  lidl  —  2udl  /m^m/^m) 
Udl  —  ^se  +  f  m(^rn/-Dm) 


(46) 


^=(1  +<*>)£?  (47) 

The  model  was  used  to  analyse  the  effect  of  the  compression 
ratio  and  the  driving  pressure  ratio  on  the  critical  entrainment 
ratio  and  to  predict  the  performance  of  an  ejector  refrigeration 
system  using  various  refrigerants.  They  found  that  the  entrainment 
ratio  is  reversely  proportional  to  the  compression  ratio,  propor¬ 
tional  with  the  driving  pressure  ratio,  and  the  system  with  R134a 
gives  better  performance. 

By  using  a  model  similar  to  that  of  Rogdakis  and  Alexis  [13],  Yu 
et  al.  [15]  derived  a  thermodynamic  model  based  on  isentropic 
process  to  estimate  the  entrainment  ratio  of  the  ejector  for  both  a 
refrigerator  and  a  pump: 


—  \  PnPmVM^py  —  /lpi )/ (fimd,s  —  ^mm))  —  1 


(48) 
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The  enthalpy  value  obtained  from  Eq.  (52)  is  corrected  for 
friction  as: 

h[  =  h[_ \  +  —  fy-l)  (57) 

Together  with  the  defining  of  Mach  number  and  sonic  velocity, 
and  given  the  initial  condition,  the  model  can  be  calculated  step  by 
step. 

The  models  were  developed  for  qualitative  and  quantitative 
studies  on  the  operation  and  performance  for  optimal  and  off 
design  operation.  The  model  can  also  be  used  to  represent  the 
global  characteristics  and  detailed  evolutions  of  all  the  relevant 
parameters.  In  addition,  it  can  be  used  to  study  the  influence  of 
individual  parameters  on  the  general  ejector  operation.  They  used 
the  model  to  analyze  the  effects  of  geometric  configurations  on 
ejector  operation  under  different  conditions  showed  that  increas¬ 
ing  ejector  diameter  improved  entrainment  ratios  and  capacities 
but  exit  pressures  decreased  accordingly.  Component  length 
affects  not  only  the  overall  ejector  size  but  also  the  capacity  and 
the  stream  mixing  process.  The  mixing  chamber  length  was 
particularly  important  for  the  control  of  the  shock  wave  intensity. 
Less  irreversibility  and  limited  superheat  was  produced  in  the 
ejector  when  the  length  was  adjusted  to  bring  supersonic  mixed 
flow  to  near  sonic  conditions  for  maximum  exit  pressure.  Hence, 
maximum  isentropic  compression  was  obtained  in  the  supersonic 
convergent,  greatly  reducing  the  shock  contribution,  such  that 
when  combined  with  the  pressure  lift  in  the  diffuser,  the  pressure 
at  the  ejector  exit  was  globally  higher. 

3.2.2.  Constant-area  mixing  model 

Keenan  and  co-workers  investigated  the  constant-area  mixing 
ejector  when  they  worked  on  the  constant-pressure  mixing  model. 
They  found  that  the  constant-area  mixing  model  offers  better 
agreement  with  the  experimental  results  than  that  of  the  constant- 
pressure  mixing  model.  Also,  the  constant-area  mixing  model  gives 
more  information  of  the  flow  inside  the  ejector  than  that  of  the 
constant-pressure  mixing  model. 

In  the  constant-area  mixing  model,  except  for  the  mixing 
method  and  the  region,  the  assumptions  made  for  the  analysis  of 
flow  within  the  ejector  are  the  same  as  those  made  in  the  case  of 
one-dimensional  constant-pressure  ejector  flow  model. 

The  models  proposed  by  Eames  et  al.  [9]  and  Huang  et  al.  [12] 
mentioned  in  Section  3.1.1  considered  that  the  constant-pressure 
mixing  occurred  in  the  constant-area  section.  Therefore,  both 
models  are  based  on  both  constant-pressure  mixing  and  constant- 
area  mixing. 

Grazzini  and  Mariani  [17]  derived  a  theoretical  model  based  on 
the  constant-area  mixing  theory.  In  their  model,  the  following 
hypotheses  were  taken: 


(1)  Ideal  gas  behavior  of  the  superheated  steam. 

(2)  Same  static  pressure  for  primary  and  secondary  flow  at  the 
nozzle  exit  plane  section. 

(3)  Constant-area,  adiabatic  mixing  process. 

(4)  Isentropic,  one-dimensional,  steady  flow  for  both  separated 
streams  and  for  the  mixed  stream. 

(5)  Losses  are  introduced  by  means  of  the  isentropic  efficiency  for 
the  primary  nozzle  and  the  diffuser. 

Using  Eq.  (16),  the  area  ratio  based  on  gas  dynamic  equation  is: 


A^ 


A 


m 


(y+i)/(y-i) 


(58) 


With  a  recursive  method  to  calculate  the  mass,  momentum  and 
energy  balance  equations  using  Mach  number  at  the  end  section, 
the  end  thermodynamic  variables  at  the  ejector  exit  can  be 
obtained  using  Eq.  (7). 

This  model  was  developed  for  designing  a  jet  pump  device  for  a 
water  refrigeration  cycle.  In  order  to  improve  the  overall 
performance  and  the  compactness  of  the  system,  they  used  the 
model  to  analyze  a  two-stage  ejector  which  showed  that  the 
solution  provides  greater  compression  ratio  while  using  the  same 
entrainment  ratio  of  the  single  stage  jet  pump.  They  also  analyzed 
a  jet  pump  with  a  third  stage  and  found  that  the  ideal  gas  model 
adopted  could  lose  physical  meaning. 

In  order  to  give  a  more  accurate  prediction  of  the  performance, 
Yapici  and  Ersoy  [18]  derived  a  local  model  based  on  constant-area 
mixing  ejector  consisting  of  a  primary  nozzle,  a  mixing  chamber  in 
cylindrical  structure  and  a  diffuser  (Fig.  6).  The  mathematical 
description  of  the  model  is  as  follows: 

Primary  nozzle:  By  using  equations,  such  as  (10),  (13),  (15), 
(49),  and  the  mass  flow  rate  at  throat  (Eq.  (1)),  they  derived  the 
following  expressions  to  describe  the  primary  nozzle  pressure 
ratio  (Eq.  (60))  and  the  nozzle  area  ratio  (Eq.  (61)) 


m 


pAu  =  PAM 


(59) 


V 

y- 1 
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1 

i  y-'1 

7lln  +  Vn{y-P 
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~MP, 
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j+VMpJ 
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Mixing  section:  For  clarity,  the  control  volume  selected  to 
analyze  the  flow  in  the  mixing  chamber  is  shown  in  Fig.  6.  From 
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Fig.  6.  Ejector  model  applied  by  Yapici  et  al. 


S.  He  et  al. / Renewable  and  Sustainable  Energy  Reviews  13  (2009)  1760-1780 


1769 


Control  surface 


Control  volume 


Fig.  7.  Control  volume  selected  to  analyze  the  flow  in  the  mixing  chamber. 


isentropic  gas  dynamic  equation,  the  static  pressure  ratio  Psl/Ppl  is, 


Psi_=  /TsoNl/ApA  /2(y,Mpi) 

Ppi  VW  VW  fi(Y, Msi) 


(62) 


where  /2  ( y,  M)  is  the  mass  flow  function: 


f2(y,M)=M 


(63) 


Based  on  the  continuity  equation  and  energy  balance  for  the 
selected  control  volume  shown  in  Fig.  7,  they  derived: 


Pm3  _  (Tm0\2(Apl\  /2(y,Mpi) 

Pp,  \TpoJ  VW  f2(Y,Mm3){  ’ 


(64) 
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PsO  =Psl  /l(y,Mpi,??n) 
PpO  PsO  /3(K>Msl) 

where 


(66) 


(67) 

Considering  that  p  =  P/(PT)  and  u2  =  M2yRT,  they  derived  the 
momentum  equation  in  the  x  direction  for  the  control  volume 
(Eq.  (68))  and  the  exit  Mach  number  of  the  mixing  chamber 
(Eq.  (69)). 
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where 


a  =  1  f4(y,Mpi)  +cof4(y,Msi) 
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Fig.  8  shows  the  control  volume  for  the  initial  interaction  region 
of  the  mixing  chamber.  This  control  volume  was  used  for  analysis 
of  the  flow  in  the  supersonic  regime.  In  order  to  analyze  the  flow  in 
this  region,  the  following  additional  assumptions  were  made: 


(1)  Streams  do  not  mix  and  are  isentropic  between  sections  1-1 
and  2-2; 

(2)  The  secondary  flow  Mach  number  at  section  2-2  is  unity, 
Ms2  =  1; 

(3)  The  primary  static  pressure  at  the  inlet  is  greater  than  that  of 
the  secondary,  Pp i  >  Ps i 

Considering  that  Ap2  =  Am 3  -Apl  and  As i  =Am 3  -Apl,  together 
with  Eq.  (61),  they  derived 
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Thus,  Mpi  can  be  found  by  using  the  above  two  equations. 
Considering  that  Ms2  =  1  and  As i  =  An3  -  Api,  they  derived  the 
inlet  static  pressure  ratio: 


(Pp2/Pp0)  (^p2/^p*  ) 

(Fpi/Fpo)  (AuAV) 


l~(Ail/An3) 

^pl/-^m3 


(i  +  yJ^i) 


Ps2  /PsO  1+7 
Psl  /FsO  A>1  /-^s2* 


(74) 


Diffuser  section:  By  applying  the  Eq.  (22)  to  calculate  the 
diffuser  pressure  ratio,  they  derived  the  following  expression  for 
the  pressure  ratio  through  the  ejector 


P40  _  P40  Pill 3  Ppl  PpO  /yrx 

PsO  Pm3  Ppl  PpO  PsO 

The  authors  derived  the  model  for  ejector  optimization. 
Meanwhile,  they  carried  out  the  parametric  study  by  applying 
the  model  and  derived  the  optimum  design  curves  for  constant- 
area  ejector  design.  Further,  they  compared  the  analytical  results 
with  similar  model  developed  by  Sun  and  Eames  [19]  which  was 
based  on  the  constant-pressure  mixing  model,  for  the  same 
operating  temperatures  of  the  ejector  refrigeration  system,  and 
found  the  optimum  coefficient  of  the  performance  and  area  ratio 
are  greater. 


3.2.  Two-phase  flow  model 


It  should  be  pointed  out  that  the  above  mentioned  models  are 
based  on  the  assumption  that  the  primary  gas  is  superheated  and 
the  flow  in  the  ejector  is  compressible  single-phase  flow  and 
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Mp2 

Pp2 


a  normal  shock  wave  may  develop  inside  the  diffuser.  However, 
under  many  application  conditions,  the  fluid  mixture  easily 
undergoes  phase  change  and  a  condensation  shock  may  develop. 
And  because  of  the  overall  rise  in  pressure  in  the  diffuser  as  well  as 
the  additional  rise  in  pressure  caused  by  the  shock  wave,  the  fluid 
mixture  is  easy  to  condensate.  On  the  other  hand,  the  secondary 
fluid  will  accelerate  to  sonic  speed  at  the  “hypothetical  throat” 
which  will  reduce  the  pressure.  Thus  the  sub  cooled  or  saturated 
secondary  fluid  will  flash  into  vapor  and  form  either  a  two-phase 
mixture  or  pure  vapor  in  the  region  occupied  by  the  secondary 
fluid  [20].  In  addition,  when  the  ejector  works  as  a  pump  using 
high  pressure  liquid  to  entrain  gas  or  higher  pressure  vapor  to 
entrain  liquid,  under  these  cases,  the  flow  is  two-phase  flow. 
Taking  account  of  these  cases,  the  models  reviewed  in  this  part  are 
capable  of  properly  simulating  ejectors  containing  two-phase 
flow. 

Sherif  et  al.  [20]  derived  an  isentropic  homogeneous  expan¬ 
sion/compression  model  to  account  for  phase  change  due  to 
expansion,  compression  and  mixing  while  attempting  to  provide 
a  more  accurate  description  of  the  constant-pressure  mixing 
process.  In  this  model,  the  primary  fluid  is  a  two-phase  mixture 
and  the  secondary  fluid  is  either  a  sub-cooled  or  saturated  liquid 
having  the  same  chemical  composition  as  the  primary  fluid. 
Therefore,  comparing  with  single-phase  flow  model,  the  biggest 
difference  in  the  model  is  the  introduction  of  dryness  of  the  fluid 
in  the  calculation  of  the  specific  volume,  enthalpy  and  entropy. 
For  clarity,  the  ejector  model  used  by  Sherif  et  al.  is  shown  in 
Fig.  9. 

Primary  nozzle:  Different  from  single-phase  flow  model,  in  the 
model,  the  primary  nozzle  is  segmented  into  two  parts,  i.e.  the 
converging  section  and  the  diverging  section.  For  given  inlet 
pressure,  temperature  and  quality,  the  specific  volume,  enthalpy 


and  entropy  can  be  calculated  by  Eqs.  (76)-(78). 


b’ni  =  E’ni.f  +Xni(l7ni)g  —  t7njf) 

(76) 

^ni  =  ^ni,f  +^ni(^ni,g  —  ^ni,f) 

(77) 

^ni  =  Snj  ,f  +^ni(Sni,g  —  Sni,f) 

(78) 

Considering  the  continuity  and 

energy  equation  in  the 

converging  section  along  with  the  definitions  of  the  Mach  number 
and  the  speed  of  sound,  the  isentropic  flow  requirement,  and  the 
equations  of  state,  the  enthalpy  and  density  are, 


Pni^ni^ni  —  Prif^nt^nt 
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By  assuming  the  pressure  at  the  nozzle  throat  Pnt  and  using 
Eq.  (81)  to  compute  the  quality  at  the  throat  xnt.  The  quality  can 
then  be  used  to  compute  the  specific  enthalpy  and  density  at  the 
throat  using  Eqs.  (82)  and  (83),  respectively.  And  then  the  Eqs.  (79) 
and  (80)  can  be  solved  in  the  velocities  um  and  unt.  The  speed  of 
sound  at  the  throat  can  be  computed  from  Eq.  (6)  and  be  compared 
to  the  velocity  term  unt.  If  the  two  quantities  are  found  different  in 
value,  another  value  of  pressure  is  tried  until  convergence.  Then 


Converging-diverging  \  Constant  pressure  \  Diffuser 
primary  nozzle  Nozzle  throat  mixing  chamber  Ejector  throat 


Fig.  9.  Ejector  model  used  by  Sherif  et  al. 
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the  mass  flow  rate  of  the  primary  fluid  can  be  calculated  by 
Eq.  (59). 

Using  the  same  method  by  assuming  the  pressure  at  the  nozzle 
exit,  the  state  parameters  at  the  nozzle  exit  can  be  solved. 

Secondary  flow:  If  the  secondary  flow  is  all  liquid,  Bernoulli’s 
equation  can  be  applied  between  the  inlet  and  exit  of  the 
secondary  flow  tube, 


,  1  2 

se  ~E  2  Pso^se 


(84) 


Assuming  that  the  liquid  density  is  nearly  constant  throughout 
the  secondary  fluid  path,  i.e.  pse  «  pso,  the  secondary  fluid  flow  rate 
can  be  calculated  using  Eq.  (59). 

For  given  entrainment  ratio,  Eqs.  (84)  and  (59)  can  thus  be 
solved  to  produce  the  pressure  and  velocity  at  the  secondary  flow 
path  exit. 

Mixing  chamber:  With  the  assumption  of  constant-pressure 
mixing,  the  continuity,  momentum,  and  energy  equations  in  the 
mixing  chamber  are 


TElp  ( 1  +  CO)  —  Pme^nie^me 


(85) 


(P  ne  —  flse)Ane  —  7Hp(l  +  CO)Ume  —  lTlpUne  TTlpOJUse 


(86) 


(hne  "E  2  Une)  ~E  (^se  +  ^  Use)^  —  (^me  +  2  Ume)0  +  &0 


(87) 


Similarly,  the  state  parameters  at  the  exit  of  the  mixing 
chamber  can  be  obtained  by  using  the  state  equations  similarly  to 
Eqs.  (81-83)  as  well  as  the  equations  defining  the  Mach  number 
and  the  speed  of  sound. 

Diffuser  section:  The  conditions  at  the  mixing  chamber  exit 
were  assumed  identical  to  those  at  the  diffuser  inlet  and  a  shock 
wave  is  assumed  to  occur,  then  the  mass,  momentum  and  energy 
balance  in  the  diffuser  are, 
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P(P,h) 

(91) 

Again,  by  assuming  the  pressure,  the  state  parameters  can  be 
calculated  using  above  equations. 

For  given  entrainment  ratio  and  inlet  fluid  states,  the  model  can 
be  used  to  determine  the  distribution  of  thermodynamic  states  as 
well  as  the  area  of  the  mixing  chamber  exit.  By  analyzing  the 
performance  of  the  ejector  applying  the  model,  the  authors  found 
that  the  variation  of  compression  ratio  with  entrainment  ratio 
follows  a  monotonically  decreasing  trend.  The  more  energetic 
primary  fluid  appears  to  be  more  suitable  for  inducing  high 
compression  ratios  than  high  entrainment  ratios. 

Cizungu  et  al.  [21]  derived  a  two-phase  thermodynamic  model 
to  calculate  the  entrainment  ratio  shown  in  Eq.  (92). 
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mentioned  above  for  simple  calculation: 
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where  y=rjexplricom  represents  the  quality  grades  of  the  jet 
compressor,  it  ranges  between  0.689  and  0.81. 

This  model  can  be  used  both  for  single-phase  flowing  ejector 
and  two-phase  flowing  ejector  with  single-component  or  two- 
component  as  the  working  fluid.  For  given  boundary  conditions  the 
optimum  geometry  of  the  ejector  can  be  calculated.  For  given 
geometry  of  the  ejector,  the  optimum  thermal  conditions  can  be 
determined.  They  used  the  model  to  analysis  the  interaction 
between  the  entrainment  ratio  and  all  relevant  dimensional 
parameters  and  compared  the  results  with  experimental  data. 
Good  agreement  was  found.  They  also  found  that  the  dimensions  of 
ejector  configuration  have  a  dominant  influence  in  deciding  the 
operating  range. 

Another  kind  of  widely  used  two-phase  flowing  ejector  is  that 
employs  liquid  as  primary  fluid  and  gas  as  secondary  fluid  or  vice 
versa.  This  kind  of  ejector  found  widely  application  in  industry, 
such  as  steam-driven  jet  pump  used  in  modern  thermal  power 
plant  to  remove  non-condensable  gases  from  the  condenser, 
ejectors  used  to  entrain  and  pump  fumes  and  dust-laden  gases 
which  otherwise  are  difficult  to  handle  and  so  on.  Therefore, 
Beithou  and  Aybar  [22]  employed  a  one-dimensional  control 
volume  method  to  develop  a  local  mathematical  model  for  the 
performance  prediction  of  the  steam-driven  jet  pump.  In  the 
model,  the  conservation  of  mass  and  momentum  equation  as  well 
as  the  entropy  equation  were  applied.  As  the  fluid  is  uncompres- 
sible,  therefore,  the  Bernoulli  equation  can  be  used  to  calculate  the 
entrained  water  velocity  at  the  nozzle.  Regarding  to  the  mixing 
section,  the  steady-state  energy  equation  was  used  to  calculate  the 
velocity  of  the  mixture  with  the  assumption  of  adiabatic  mixing 
process  as  well  of  no  potential  energy  change.  Take  the  head  loss 
into  account,  the  Bernoulli  equation  was  used  in  the  diffuser  to  find 
out  the  pressure  and  velocity  profiles  along  it.  The  mathematical 
model  is  shown  as  follows. 

Steam  nozzle: 

Conservation  of  mass  :  ^PjUjAj  =  ^peiieAe  (1) 


Conservation  of  momentum 


du  du  dP  _ 
P  ~Fr~  +  p  tl  — — b  — —  0 


dt 


dx  dx 


(94) 


The  governing  equations  of  steady-state  subsonic  or  supersonic 
isentropic  flow: 


dx 


(j puA )  =  0 


(95) 


du  1  dP 

U  dx  +  p  dx 


=  0 


(96) 


Also  the  entropy  equation  Eq.  (11)  is  used. 

Water  nozzle: 

By  using  Bernoulli  equation  between  the  water  tank  and  the 
exit  of  the  water  nozzle  and  the  continuity  equation,  the  water 
velocity  and  water  flow  rate  can  be  calculated  by  Eqs.  (97)  and 
(59): 


(97) 


where  fm  represents  the  friction  factor,  and  £m  is  the  momentum 
loss  factor,  the  enthalpy  differences  are  given  by  A/iexp  =  hp0  -  hp i>is 
pi  is  and  A/iexs  =  hs>t  -  hsi  iS.The  authors  simplified  the  equations 


where,  the  subscript  “T”  defines  the  water  properties  in  the  supply 
water  tank,  and  the  subscript  “wc”  shows  the  water  properties  at 
the  inlet  of  water  nozzle. 
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Mixing  section:  The  steady-state  energy  equation  with 
assumption  of  no  potential  energy  change  is  similarly  to  Eq.  (30). 

Diffuser  section:  Applying  Bernoulli  and  considering  the  head 
loss,  they  derived: 


i!d  +  Pd 
2  Pd 


+  hi 


where  hL  is  the  head  loss  which  is  given  by, 
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where  rjp  is  the  pressure  recovery  coefficient. 

In  the  steam  nozzle,  the  iterative  finite  difference  procedure 
with  under-relaxation  is  used  to  approximate  the  derivatives  in  the 
steam  nozzle  governing  equations  that  are  expressed  in  terms  of 
the  mesh  points. 

The  initial  condition  is  the  state  parameters  and  geometrical 
parameter  of  the  inlet  section  of  the  steam  subsonic  nozzle.  The 
nozzle  shape  is  such  that  the  area  variation  with  distance  along  the 
nozzle  is  given  by  a  quadratic  function.  In  the  subsonic  nozzle,  the 
effect  of  the  density  changes  is  negligible.  Assuming  p,  =  pA  gives 
the  initialization  of  density  at  every  grid  point,  and  then  calculates 
the  velocity  and  the  pressure  at  each  grid  point  for  the  first 
iteration  by, 
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where  the  index  1  shows  the  initial  condition,  n  shows  the  density 
at  the  previous  iteration,  and  the  under-relaxation  factor  r  is  to  be 
chosen  less  than  1.  Then  the  pressure,  velocity  and  density  can  be 
updated  by  following  equations: 


Pi 


(103) 


Fig.  10.  Control  volume  for  the  modeling  of  the  mixing  chamber. 


consideration,  Deberne  et  al.  [24]  derived  a  global  volume  model 
(Fig.  10)  for  steam  ejector  performance  with  the  consideration  of 
void  fraction,  area  contraction  ratio  and  condensation  rate  in  the 
conservation  of  mass,  momentum  and  energy  equations  based  on 
steady  and  one-dimensional  assumption. 

By  defining  the  following  parameters: 


Void  fraction  :  = 

Area  contraction  ratio  : 


1 


P  lsuls 

P\\u\\ 


Q 


As  +  Ai 


(106) 

(107) 


Condensation  rate  :  R  = 

The  conservation  equations  applied  to  the  global  volume  lead 
to: 


X\LV  \ 

CpiAr^y 


(108) 


Mass:  p2u2  =  (1  +  co)  -^SiPisiq. 


(109) 


Momentum  :  (1  +  co)u2  + 


Q 


—  (uls  +  OJU^is)  + 


SlsPisAs 

Pis 


P2+  1 


h) 


£\  PlsUls 


(110) 


Energy  :  (1  +  co) 

'  u22] 

=  CO 

ii2 

h  J_  11 

hu+  2 

+ 

u  ,<} 
fi,s+  2 

State  equation  : 

h2  =  h{P 

2 1P2) 

(ill) 

(12) 


U  i  =  U[_ -[ 


-(Pi- Pi-1 )- 

Pi  Lli 


(104) 


p"  +  r 


Pi  A  A 
U[A[ 


(105) 


The  goal  of  the  model  is  to  evaluate  the  pressure,  density  and 
velocity  distribution  along  the  ejector.  The  heat  transfer  rate  from 
steam  to  water  was  obtained  using  the  direct-contact  heat  transfer 
coefficient  as  a  function  of  the  length  of  a  conical  steam  jet  plume. 
The  condensation  profile,  non-condensed  steam  percentage  in 
volume  versus  non-dimensional  length  of  mixing  section  was 
obtained  from  the  heat  transfer  rate  profile.  They  compared  the 
calculated  pressure  profile  of  the  model  with  Cattadori  et  al.’s  [23] 
experimental  pressure  profile  which  found  that  the  experimental 
and  calculated  pressure  distributions  are  in  good  agreement  in  the 
mixing  nozzle  and  diffuser,  and  there  are  some  difference  in  steam 
nozzle. 

The  above  mentioned  local  model  models  the  flow  in  the 
injector  at  each  point,  and  taking  into  account  major  phenomena 
such  as  condensation  of  droplets,  heat  and  momentum  transfers, 
viscous  dissipation  on  the  wall  non-adiabatic  and  non-equilibrium 
flow.  Such  a  model  is  more  realistic,  but  needs  a  perfect  knowledge 
of  the  previously  mentioned  phenomena  as  well  as  their  mutual 
interactions  that  makes  the  model  very  complex.  Take  these  in 


where  the  equivalent  pressure  P*  are  obtained  by  the  empirical 
correction  numerically  introduced  from  the  experiment  results  at 
the  CETHIL: 

P*  =  1.01R2  572  x  Pls  (112) 

wherein,  the  subscript  of  Is  refers  to  mixing  section  inlet.  The 
subscripts  11  and  2,  respectively,  denotes  the  water  mixing  section 
inlet  (water  nozzle  outlet),  the  mixing  section  outlet  planes. 
AT  =  Tsat  -  Tu  (Tsat  is  saturation  temperature  at  pressure  P)  is  the 
inlet  liquid  sub  cooling,  Lv  is  the  latent  heat  of  water,  Cpi  is  the  heat 
capacity  of  liquid  and  is  the  vapor  quality  at  the  mixing  inlet 
section. 

Shock  wave: 

They  discussed  the  shock  wave  under  two  conditions:  (i)  too 
high  void  fraction,  the  fluid  will  not  be  condensed  enough  to  reach 
a  complete  condensation  across  the  shock  wave  and  (ii)  the  void 
fraction  is  too  low,  the  fluid  velocity  will  be  too  low  (lower  than  the 
sound  speed)  and  no  shock  wave  can  occur. 

The  first  sound  speed  corresponds  to  the  case  where  response  of 
droplets  is  negligible  (droplets  are  frozen  with  no  mass  and 
momentum  transfers),  and  is  similar  to  the  sound  speed  in  pure 
gases: 

Cf  =  \J  yRgTs  (113) 


S.  He  et  al / Renewable  and  Sustainable  Energy  Reviews  13  (2009)  1760-1780 


1773 


The  other  sound  speed,  presented  in  Young  and  Guha  [25]  is 
called  the  full  equilibrium  sound  speed,  where  all  processes  are  in 
equilibrium  during  shock  wave: 


ce  = 


xyRgTs 


(114) 


where  c  =  cps  +  ((1  -  e)/e)cp i,  in  the  model,  the  normal  shock  wave 
is  consider  in  full  normal  shock  wave,  fully  dispersed  with 
discontinuity  in  flow  properties  and  partly  dispersed  shock  wave. 
For  shock  wave,  and  neglect  the  effect  of  the  wall  on  the  wave, 
surface  area  of  sections  2-2  and  3-3  are  taken  to  be  equal,  and  then 
the  conservation  equations  are  as  follows: 


Mass:  p3u3  =  p2u2  with  p2  =  sp2s  +  (1  -  s)p2\  (115) 


Momentum  :  p3uj  +  P3  =  p2u2  +  P2 


(116) 


Energy:  p3u3  (h3  +  lufj 

uP 

=  /o2^-  +  u2[e2P2s,I2s  +  (l  —  e2)y02ift2i]  (117) 

State  equation  :  h3  =  h(P3,  p3)  (12) 

Diffuser:  By  introducing  a  pressure  loss  term  in  the  Bernoulli 
equation,  then  the  Bernoulli  equation  is, 

P4+2P3u24  =  P3+^p3uj-ridp3uj  (118) 

Once  the  inlet  conditions  are  known,  as  well  as  the  steam 
injector  geometry,  the  condensation  rate  R  is  calculated.  The  term 
P*  is  calculated  with  the  equivalent  pressure  ratio  r*  and  then 
all  the  other  constants.  If  no  solution  for  the  system  of  Eqs. 
(116)  +  (12)  +  (117)  +  (118)  can  be  found,  the  calculation  aborts  if 
the  flow  velocity  u2  is  lower  than  sound  speed  upstream  of  the 
shock  wave  to  ensure  the  assumption  (two-phase  supersonic 
upstream  flow,  shock  wave  and  liquid  downstream  flow).  Physical 
quantities  of  the  flow  are  then  calculated  upstream  and  down¬ 
stream  of  the  shock  wave  and  the  minimum  void  fraction  condition 
is  a  posteriori  verified  (to  be  sure  that  a  normal  shock  wave  exists). 
The  outlet  quantities  are,  in  this  case,  calculated. 


The  authors  also  set  up  a  test  rig  to  verify  the  model  and  15% 
performance  prediction  accuracy  had  been  found.  Then  a  series  of 
parametric  investigation  were  carried  out  by  applying  the  model  to 
predict  the  influence  of  different  geometrical  parameters  of  the 
ejector  (mixing  section  throat  diameter,  position  of  the  liquid  pipe, 
etc.)  and  the  influence  of  physical  parameters  (temperature,  liquid 
and  steam  pressure,  etc.)  on  the  ejector  performance. 

The  thermodynamic  models  summarized  in  this  part  are 
reviewed  in  Table  1. 

4.  Dynamic  model 

Despite  the  usefulness  and  the  remarkable  progress  the 
thermodynamic  models  provided  for  the  general  understanding 
of  ejectors,  this  kind  of  models  were  unable  to  correctly  reproduce 
the  flow  physics  locally  along  the  ejector.  It  is  the  understanding  of 
local  interactions  between  shock  waves  and  boundary  layers,  their 
influence  on  mixing  and  recompression  rate,  that  will  allow  a  more 
reliable  and  accurate  design,  in  terms  of  geometry,  refrigerant  type 
and  operation  conditions.  A  way  of  achieving  this  objective  at  a 
reasonable  cost  is  through  computational  fluid  dynamics  (CFD) 
modeling  approach  which  provides  a  better  understanding  of  the 
hydrodynamics  of  the  ejector.  Such  models  are  called  dynamic 
model  (listed  in  Table  2).  Although  dynamic  model  needs  more 
time  compared  with  thermodynamic  model,  it  is  more  related  to 
the  true  physics  of  the  problem.  Therefore,  dynamic  model  is  likely 
to  give  better  agreement  with  experiment  results.  With  the 
development  of  computer  science  and  mathematical  method, 
more  and  more  researchers  have  focused  on  this  method  with 
attempt  to  obtain  more  accurate  information  in  the  ejector. 

4.2.  Single-phase  flow  model 

The  single-phase  flow  model  includes  gas-gas  ejector  and 
liquid-liquid  ejector  flow  model.  The  former  considers  the  flow  as 
compressible  while  the  latter  as  incompressible. 

Because  of  poor  mesh  resolution,  the  early  CFD  models  [26-28] 
are  failed  to  track  the  shock  wave  and  poor  agreement  with  the 
experimental  results.  Rusly  et  al.  [29]  simulated  the  flow  through 
an  R141b  ejector  by  using  the  real  gas  model  in  the  commercial 
code,  FLUENT.  The  effects  of  ejector  geometries  were  investigated 


Table  1 

Features  of  thermodynamic  model  reviewed  in  present  work. 


Ref. 

Media 

Gas  property 

Ideal  gas 

Multiphase 
flow  model 

Real  gas 

Mixing  section 

Mixing  mechanism 

NXPa  Section  y-y 

CPMb 

CAMC 

B.C.d 

I.C.e 

Other  relations 

[9] 

Steam 

V 

Compressible  single-phase 

si 

V 

V 

Tp,  Pp,  mp,  Ts,  Ps,  ms 

Pl/Ps.t 

Eqs.  (15-21) 

[10] 

Steam 

V 

Compressible  single-phase 

V 

Tp,  Pp,  Ts,  Ps 

CO 

Eqs.  (15),  (23-28) 

[12] 

R141b 

V 

Compressible  single-phase 

V 

V 

V 

Tp,Pp,Ts,Ps,P*,At 

Ant>  Ai 

Eqs.  (6-10),  (13),  (29),  (30) 

[13] 

R717 

V 

Compressible  single-phase 

V 

V 

r,,  t2,  r3 

Eqs.  (31-41) 

[14] 

R134a* 

V 

Compressible  single-phase 

V 

V 

Tp,  Pp,  Ts,  Ps,  Pe 

Eqs.  (42-47) 

[15] 

R142b 

V 

Compressible  single-phase 

V 

V 

Inlet  conditions 

Eqs.  (48-51) 

[16] 

R142b 

V 

Compressible  single-phase 

V 

V 

Inlet  conditions 

p 

Eqs.  (5),  (6),  (52-57) 

[17] 

Steam 

V 

Compressible  single-phase 

7 

V 

V 

Tp,  Pp,  Ts,  Ps 

Eqs.  (l)-(3),  (7),  (16),  (58) 

[18] 

R123 

7 

Compressible  single-phase 

V 

V 

V 

Tp,  Ts,  Te 

Eqs.  (10),  (13),  (15),  (49), 

(59-75) 

[20] 

R134a 

V 

Homogeneous  model 

V 

Tp,  Pp,Ts,  Ps 

Pnt,  Pyl 

Eqs.  (6),  (59),  (76-91) 

[21] 

nh3 

V 

Two-phase 

Tp,  Ts,  Te 

Eqs.  (92),  (93) 

[22] 

Steam-water 

Condensation 

Inlet  conditions 

Eqs.  (1),  (11),  (30),  (94-105) 

[24] 

Steam 

Condensation 

V 

Inlet  conditions 

Eqs.  (12),  (109-118) 

a  NXP:  nozzle  exit  plane. 
b  CPM:  constant-pressure  mixing. 
c  CAM:  constant-area  mixing. 
d  B.C.:  boundary  conditions. 
e  I.C.:  initial  conditions. 

Other  media  used  are  R152a,  R290,  R600a,  R717. 
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Table  2 

Features  of  dynamic  model  reviewed  in  present  work. 


Ref. 

Media 

Gas  property 

Multiphase  flow  model 

Turbulence  model 

Mixing  mechanism 

Other  relations 

Ideal  gas  Real  gas 

CPM 

CAM 

B.C. 

[29] 

R141b 

V 

Compressible  single-phase 

Realizable  k-e 

V 

Tp,  Ts,  Te 

Eqs.  (4),  (119-122) 

[31] 

Air 

V 

Compressible  single-phase 

RNG  k-e,  SST,  /<-(*) 

V 

Tp,  Pp,  Ts,  Ps,  Pe 

Eqs.  (4),  (119-126) 

[34] 

R142b 

V 

Compressible  single-phase 

k- e,  SST,  k- co 

V 

Tp,  Pp,  Ts,  Ps,  Pe 

Eqs.  (4),  (119-126) 

[35] 

R141b 

V 

Compressible  single-phase 

RNG  k-e 

V 

Pp,  Ps,  Pe 

Eqs.  (4),  (119-124) 

[36] 

Steam 

V 

Compressible  single-phase 

Realizable  k-e 

V 

V 

Pp,  Ps,  Pe 

Eqs.  (4),  (119-124) 

[37] 

Water-air 

Air 

ASM 

Standard  k-e 

Up,  Ps,  Pe 

Eqs.  (124),  (127-130) 

[38] 

Water-air 

Air 

Mixture 

Up,  Ps,  Pe 

Eqs.  (127-130) 

[39] 

Water-air 

Air 

ASM  Eulerian 

Standard  k-e 

Up,  Ps,  Pe 

Eqs.  (127-133) 

[41] 

Air-water 

Air 

mixture 

Standard  k-e 

u?p,  ms,  Pe 

Eqs.  (127-130) 

[42] 

Water-vapor 

Separated  two-phase 

Pp,  Tp,  Ps,  x 

Eqs.  (134-149) 

numerically.  The  CFD’s  results  were  validated  with  experimental 
data  provided  by  others  [12]  and  good  agreements  were  found. 

As  it  is  often  the  case  in  transonic  compressible  flows  involving 
shock  reflections  and  shock-mixing  layer  interaction,  the  choice  of 
the  turbulence  model  and  grid  refinements  are  critical  points  [30]. 
Therefore,  Bartosiewicz  et  al.  [31]  used  the  CFD  model  aiming  at 
validating  the  choice  of  a  turbulence  model  for  the  computation  of 
supersonic  ejectors  in  refrigeration  applications.  In  order  to  reduce 
the  complexity  of  the  model  and  to  be  able  to  use  the  available 
experimental  data,  they  chose  air  as  the  working  fluid.  Six 
turbulence  models,  named  k- s,  realizable  k- s,  renormalization- 
group  (RNG)  k- 8,  RSM  and  shear-stress  transport  (SST)  k- to  have 
been  tested  and  compared  with  measurements  of  Desevaux  and 
Aeschbacher  [32]  and  Desevaux  et  al.  [33].  The  compressible 
steady-state  axisymmetric  form  of  the  fluid  flow  governing 
equations  are  shown  as  Eqs.  (119-122)  together  with  ideal  gas 
law  Eq.  (123). 


Conservation  of  mass  : 


Adp  + 
A9t  + 


0 


(119) 


Conservation  of  momentum  :  ^  ( pux )  +  (pzijZij) 

=  _9P  9ry 
dX  [  ~  dXj 


(120) 


Conservation  of  energy  :  —  ( pE )  +  —  -(Ui(pE  +  P)) 


9T  - 

—  V(Q'eff  ^  )  +  V(Zij(Tjj)) 

(121) 

with 

Ty  Meff  ( 

^9zij  9zij\ 
v  9xj  1  9xj ) 

2  9uk  „ 

3  Meff  dxk  d,J 

(122) 

P 

P~  RT 

(123) 

As  the  k- s,  RNG  k-e  and  k- to  models  are  based  on  the  Boussinesq 
hypothesis  which  makes  these  approaches  need  relatively  low 
computational  cost  associated  with  the  determination  of  the 
turbulent  viscosity.  However,  the  main  drawback  of  the  approach 
is  the  assumption  that  the  turbulence  is  isotropic.  While  the 
Reynolds  stress  model  (RSM)  does  not  rely  on  the  Boussinesq 
assumption,  but  the  associated  CPU  cost  is  relatively  high. 
Comparing  the  vast  simulations  by  applying  these  approaches, 
Bartosiewicz  and  his  coworkers  found  that  the  RNG  k-e  and  SST 


k- co  models  were  very  promising  for  ejector  analysis.  Therefore, 
they  compared  these  two  approaches  in  detail. 

The  RNG-k-s  model  (RNG) 


+  Cu^Gk-C2ep^-Re  (124) 

With  Rr  =  il  =  Sk/s,  r/o  =  4.38,  £  =  0.012, 

Ci8  =  1.42,  C2e  =  1.68,  C^  =  0.0845,  and  /xt  =  pC|JL(/<2/e). 

The  standard  and  SST  k-co  models 


9 

dX[ 


(pm) 


9 


dx\ 


OtsP'eff 


ds 


dx; 


1J 


9  9  9 

9t(p®)  +  ax[(p"Ui)  =  dxj 

With  fit  =  8*  ^  and 


[jj  + 


lJLt\  dco 

9*jJ 


o 


+  G, 


CO  1  CO 


(125) 


d  3 

gp  (PCO)  +  (pCOUj) 


-  Ya  +  (1 


(126) 


The  investigation  has  shown  that  the  RNG  k-s  and  SST  k- co 
models  were  the  best  suited  to  predict  the  shock  phase,  strength, 
and  the  mean  line  of  pressure  recovery  with  the  SST  k-co  model 
shown  better  performances  in  term  of  stream  mixing. 

Then  they  used  the  validated  model  to  reproduce  the  different 
operation  modes  of  a  supersonic  ejector,  ranging  from  on-design 
point  to  off-design  to  capture  the  shock-boundary  layer  interac¬ 
tions,  such  as  boundary  layer  separation,  flow  separation  and 
recirculation  at  the  motive  nozzle  exit  and  the  diffuser.  The  results 
turned  out  that  CFD  is  an  efficient  diagnosis  tool  of  ejector  analysis 
(mixing,  flow  separation),  for  design,  and  performance  optimiza¬ 
tion  (optimum  entrainment  and  recompression  ratios). 

In  Bartosiewicz  et  al.’s  [34]  work,  they  improved  the  previously 
mentioned  model  by  using  real  gas  thermodynamic  and  trans¬ 
portation  properties  derived  from  NIST-REFPROP  database  to 
perform  simulation  studies  on  ejector  operation  with  R142b.  They 
observed  over-expanded  and  a  strong  oblique  shock  at  the 
primary  nozzle  exit  which  induces  a  mixing  layer  separation 
associated  with  a  corresponding  energy  loss.  Consequently,  this 
separation  creates  an  adverse  pressure  gradient  at  the  secondary 
inlet,  giving  rise  to  favorable  conditions  for  an  outflow.  The  one¬ 
dimensional  model  cannot  take  into  account  such  irreversible 
losses  induced  by  oblique  shock  wave  reflection  and  flow 
separation.  This  local  CFD  modeling  takes  into  account  shock¬ 
boundary  layer  interactions  in  a  real  refrigerant.  The  numerical 
results  contribute  to  understanding  the  local  structure  of  the  flow 
and  demonstrate  the  critical  role  of  the  secondary  nozzle  for  the 
mixing  rate  performance. 
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Zhu  et  al.  [35]  employed  the  CFD  technique  with  RNG  k- s  model 
to  investigate  the  effects  of  two  important  ejector  geometry 
parameters:  the  primary  Nozzle  Exit  Position  (NXP)  and  the  mixing 
section  converging  angle  0,  on  its  performance.  The  flow  inside  the 
ejector  is  governed  by  the  compressible  steady-state  turbulent 
form  of  the  flow  governing  equations.  They  found  that  the 
optimum  NXP  is  not  only  proportional  to  the  mixing  section  throat 
diameter,  but  also  increases  as  the  primary  flow  pressure  rises.  On 
the  other  hand,  the  ejector  performance  is  very  sensitive  to  0 
especially  near  the  optimum  working  point.  The  entrainment  ratio 
varies  as  much  as  26.6%  by  changing  6.  In  order  to  better  maximize 
the  ejector  performance,  a  relatively  bigger#  is  required  when  the 
primary  flow  pressure  rises.  These  findings  are  good  for  guiding  the 
adjustment  of  NXP  and  6  in  order  to  obtain  the  best  ejector  system 
performance  when  the  operating  conditions  are  different  from  the 
on-design  conditions. 

The  previously  mentioned  dynamic  models  are  based  on 
constant-pressure  mixing  ejector.  Pianthong  [36]  employed  the 
CFD  with  realizable  k- s  turbulence  model  to  predict  the  flow 
phenomena  and  performance  of  both  constant-pressure  mixing 
and  constant-area  mixing  steam  ejectors  with  application  in 
refrigeration  system.  The  results  indicate  that  CFD  can  predict 
ejector  performance  very  well  and  reveal  the  effect  of  operating 
conditions  on  the  effective  area  that  is  directly  related  its 
performance.  Besides,  they  found  that  the  flow  pattern  does  not 
depend  much  on  the  suction  zone  because  the  results  of 
axisymmetric  and  3D  simulation  are  similar. 

4.2.  Two-phase  flow  model 

There  are  two  approaches  for  the  numerical  calculation  of 
multiphase  flows:  the  Euler- Lagrange  approach  and  the  Euler- 
Euler  approach.  In  the  latter  approach,  different  phases  are  treated 
as  interpenetrating  continua  and  this  approach  has  been  widely 
adopted  in  ejector  modeling.  The  conservation  equations  have 
similar  form  for  all  phases  in  this  approach.  There  are  two  Euler- 
Euler  multiphase  models  for  ejector  modeling:  the  mixture  model 
and  the  Eulerian  model. 


4.2.1.  Mixture  model 

This  model  solves  the  momentum,  continuity,  and  energy 
equations  for  the  mixture  while  the  volume  fraction  equations  are 
solved  only  for  the  dispersed  phases.  The  model  uses  algebraic 
expressions  for  the  relative  velocities  and  the  equations  for 
multiphase  flows  are  derived  by  mass-weighted  averaging  or 
Favre-averaging  method. 

The  continuity  and  the  momentum  equations  for  mixture  flow 
are, 


(Q'mPm)  +  V(Q'mpmUm)  —  r m 


(127) 


d 

(pmum)  +  V (pmumum)  =  —  VP  +  V(Tm  +  tT m)  +  V rDm 

+  (128) 

The  diffusion  stress  term  \yrDm  represents  the  momentum 
diffusion  due  to  the  relative  motions  of  the  two  phases  while  Mm  is 
the  influence  of  the  surface  tension  force  on  the  mixture  and 
depends  on  the  geometry  of  the  interface. 

In  order  to  calculate  the  relative  velocity,  it  is  necessary  to 
calculate  the  average  interfacial  momentum  source  for  the 
dispersed  phase  Mp.  By  neglecting  the  surface  tension  forces 
(Mm  =  0)  and  assuming  that  the  phase  pressures  are  equal  and 


rearranging,  the  equation  for  dispersed  phase  Mp  is 


Mr 


^p(Pp  Pm) 


du 


m 


dt 


+  (umV)um  —  g 


dUMp  , 

+  apPp  - E  (uMpV)Umf 

+  jo'pPp  _(UmV)UMp  +  (^Mp  V)Ump]  | 

+  {— V  _cyp(Tp  +  fjp)]  }  +  otp'V (Tm  +  zTm  +  rDm) 


(129) 


The  second  term  from  Eq.  (128)  is  omitted  by  assuming  local 
equilibrium  (the  bubbles  are  rapidly  accelerated  to  the  terminal 
rise  velocity).  The  third  term  corresponds,  in  rotational  motion,  to 
the  Coriolis  force  while  the  last  term  is  the  contributions  of  the 
viscous,  turbulent  and  diffusion  stresses.  These  are  usually 
neglected. 

By  neglecting  the  virtual  mass  effect,  lift-force  effect,  basset  (or 
history)  effect  except  the  viscous  drag  and  assuming  spherical 
bubbles,  the  momentum  transfer  induced  by  the  fluid-fluid 
interaction  force  is 


Mr 


—  pcApCo  |  ucp  |  uCp 


T/p(pp  Pm) 


-g  +  (limV)Um  + 


du 


m 


dt 


(130) 


where  ucp  is  the  relative  velocity,  ucp  =  uc-  up. 

Kandakure  et  al.  [37]  used  the  algebraic  slip  mixture  (ASM) 
model  to  study  the  hydrodynamic  characteristics  of  the  ejector. 
The  model  accounts  for  the  relative  velocity  between  the 
continuous  and  dispersed  phases  as  the  interfacial  momentum 
exchange  term.  In  the  simulation,  water  was  taken  as  the  primary 
fluid  and  air  as  the  entrained  fluid.  The  air  phase  was  assumed  to 
obey  the  ideal  gas  law.  Standard  k-&  model  was  used  for  modeling 
the  turbulent  behavior  of  the  flow.  The  slip  between  the  phases 
was  taken  as  some  percentage  of  local  primary  velocity.  They  found 
that  at  low  value  of  area  ratio  (ratio  of  the  mixing  chamber’s  throat 
area  to  nozzle  throat  area),  due  to  the  larger  diameter  of  the  water 
jet,  the  annular  area  available  for  air  flow  reduces,  causing 
recirculation  of  the  entrained  air  within  the  converging  section  of 
the  ejector.  On  the  other  hand,  for  higher  values  of  area  ratio,  due  to 
smaller  diameter  of  the  water  jet,  the  momentum  transfer  to  the 
air  decreases  and  all  the  entrained  air  cannot  be  forced  through  the 
throat.  As  a  result,  the  net  air  flow  rate  going  into  the  throat  for 
both  area  ratios  is  small.  Thus  there  is  an  optimum  area  ratio  for 
the  maximum  air  entrainment  rate.  The  air  entrainment  rate 
correlates  with  pressure  difference  between  the  air  entry  and 
throat  exit  for  a  wide  variety  of  ejector  geometries  and  operating 
conditions.  The  overall  head  loss  factor  and  the  ejector  efficiency 
can  be  predicted  a  priori. 

Kim  et  al.  [38]  also  investigated  the  hydrodynamic  character¬ 
istics  of  the  liquid-gas  (water-air)  ejectors  with  mixture  model 
and  experiments  focusing  on  the  gas  hold-up  performance.  They 
found  the  gas  phase  hold  up  increases  with  increasing  liquid 
circulating  rate  but  it  decreases  with  increasing  liquid  level  in  the 
column.  The  gas  suction  rate  increases  with  increasing  the  liquid 
circulating  rate  but  it  decreases  with  increasing  the  liquid  level  in 
the  column  and  nozzle  diameter. 

In  order  to  get  maximum  entrainment  of  the  secondary  fluid, 
Yadav  and  Patwardhan  [39]  used  the  ASM  model  to  determine  the 
optimum  geometry  of  suction  chamber  of  the  ejector  which  is 
worked  as  a  contactor  for  two  fluids  mixing.  It  is  found 
experimentally  that  the  entrainment  of  the  secondary  fluid  was 
more  when  straight  portion  of  the  throat  is  of  zero  length  [40]. 
Hence,  the  geometrical  model  of  the  ejector  adopted  by  Yadav  and 
Patwardhan  is  shown  in  Fig.  1 1  with  zero  length  of  throat.  Also, 
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water  was  taken  as  the  primary  fluid  and  air  as  the  secondary  fluid 
with  air  assumed  to  obey  the  ideal  gas  law.  The  standard  /c-s 
turbulence  model  (per  phase)  was  used  to  calculate  the  turbulence. 
Then  they  carried  out  a  large  number  of  CFD  simulations  using  the 
model  to  understand  the  effects  of  the  geometrical  parameters  of 
the  suction  chamber  on  entrainment  of  the  secondary  fluid.  They 
found  that  the  geometry  of  the  suction  chamber  has  significant 
effect  on  the  rate  of  entrainment  of  the  secondary  fluid.  After 
studying  the  three  parameters:  the  projection  ratio  (PR:  distance 
between  nozzle  tip  and  entry  to  throat),  area  ratio  and  angle  of  the 
converging  section  based  on  the  model,  Yadav  and  Patwardhan 
recommended  that  a  PR  of  five  produces  the  highest  rate  of 
entrainment  of  the  secondary  fluid.  Increasing  the  PR  beyond  five 
does  not  lead  to  any  significant  improvement  in  the  suction 
capacity  of  the  jet.  In  addition,  they  found  a  suction  chamber  with 
(PI  -  Dl)/Dl  =  6.6  (Ds:  diameter  of  suction  chamber,  Dn:  diameter 
of  nozzle  exit)  gives  the  largest  entrainment.  A  larger  suction 
chamber  leads  to  a  reduction  in  suction  capacity.  The  angle  of 
converging  section  is  recommended  in  the  range  of  5°-15°. 

The  previously  mentioned  models  are  taken  water  as  primary 
fluid  and  air  as  secondary  fluid.  Balamurugan  et  al.  [41  ]  used  the 


mixture  model  with  standard  k-e  turbulence  model  to  investigate 
the  hydrodynamic  characteristics  of  ejectors  using  air  as  the 
primary  fluid  and  water  as  the  entrained  fluid.  Both  air  and  water 
were  considered  at  room  temperature.  The  air  was  assumed  to 
obey  an  ideal  gas  law.  They  observed  that  the  liquid  entrainment 
rate  increases  with  an  increase  in  the  liquid  level  and  the  gas 
velocity.  The  entrainment  rate  was  found  to  be  the  highest 
corresponding  to  area  ratio  of  about  0.155. 

4.2.2.  Eulerian  model 

The  Eulerian  model  solves  the  momentum  and  continuity 
equations  for  each  phase  separately.  The  coupling  of  the  equations 
is  achieved  through  the  pressure  and  inter-phase  exchange 
coefficients.  Momentum  exchange  between  the  phases  is  also 
dependent  upon  the  type  of  mixture  under  consideration. 

The  continuity  equation  and  the  momentum  equations  for  the 
Eulerian  model  are  the  same  as  the  mixture  model.  The  average 
interfacial  momentum  source  for  the  phase  k  is 


—  J2(Rjk  +  rnjkUjk)  +  otkPk(Fk  T-  Eiift,/<  +  Fvm j<)  (131 ) 

j= i 


where  Fk  is  an  external  body  force,  FUfttk  is  a  lift  force,  Fvm,/<  is  a 
virtual  mass  force,  Rjk  is  an  interaction  force  between  phase  i  and 
phase  k  and  p  is  the  pressure  shared  by  all  the  phases.  The  lift  force 
acts  on  a  dispersed  phase  due  to  velocity  gradient  in  the 
continuous-phase  flow  field,  is: 

Fiift p  =  -0.5pcap(uc  -  up)  x  (V  x  uc)  (132) 

For  multiphase  flows,  the  virtual  mass  effect  occurs  when  a 
dispersed  phase  accelerates  relative  to  the  continuous  phase.  The 
inertia  of  the  continuous  phase  mass  encountered  while  accel¬ 
erating  dispersed  phase  exerts  a  “virtual  mass  force”  on  the 
dispersed  phase,  which  is  given  by: 


Fym,p  —  0.5Qfp/7c(^^  (Uc  —  Up)) 


(133) 


Eulerian  model  assumes  dispersion  type  flow  where  the 
secondary  fluid  is  distributed  in  the  primary  fluid.  However,  in 
Yadav  and  Patwardhan’s  [39]  work,  the  ejector  was  arranged  such 
that  the  water  jet  was  directed  downwards.  The  entrained  air 
flowed  around  the  water  jet  in  the  annular  space  between  the 
water  jet  and  the  ejector.  They  found  no  bubble  formation  inside 
the  ejector  and  both  the  fluids  flowed  co-currently.  Since  there  was 
no  bubble,  the  virtual  mass,  which  arises  from  acceleration,  and  lift 
force,  which  arises  from  unbalanced  pressure  forces  around  the 
bubble,  do  not  exist.  Hence,  the  Eulerian  model  is  not  a  suitable 
choice  for  studying  the  coaxial  type  of  flows  in  ejector. 

The  previously  mentioned  models  consider  the  primary  and 
secondary  fluids  are  of  different  chemical  compositions  with  non¬ 
condensable  gas.  However,  when  both  fluids  are  of  the  same 
chemical  composition,  the  vapor  is  condensable  in  the  fluid,  there 
may  be  need  a  different  model  to  simulate  the  hydrodynamic 
performance  of  the  ejector.  Therefore,  Narabayashi  et  al.  [42] 
derived  a  separate  two-phase  flow  model  to  analysis  the 
performance  of  a  steam  injector  entrained  by  high  pressure  water 
in  order  to  check  the  feasibility  of  a  large-scale  steam  injector  for 
which  a  demonstration  test  was  not  able  to  be  conducted.  In  the 
model,  attention  was  paid  on  the  radial  heat  transport  enhance¬ 
ment  by  the  turbulence  because  of  relatively  large-scale  water  jet 
being  used.  Also,  the  interface  activation  was  taken  into  account. 
The  model  including  steam  phase,  water  phase  and  the  interface  is 
shown  as  follows. 
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Steam: 


5.  Empirical/semi-empirical  model 


div(PgOg)  =  -mg  | 
div(pgL/gVg  -  TuggradVg)  = 

div(pgL?gWg  -  FuggradWg) 

div(pgL?g/ig  -  F/iggrad/ig)  = 
div(pgUgkg  -  F/Cggrad/<g)  = 

div(pgUg6g  -  rgggradgg)  = 

Water 

div(p1L?1)  =  -mgi 
div(p1L?1\/1  -  FuigradVi)  =  - 


(134) 


9P 

dr 

(135) 

3P 

=  -Q z  +  °i  Wg-Wi\(Ws 

-wo 

(136) 

DP  ^ 

+  <p 

DT 

(137) 

Gkg  —  PgSg 

(138) 

£ 

C,  jfCkg  -  C2pgsg/kg 

ICg 

(139) 

(140) 

dP 

dr 

(141) 

di v(p1L?1W1  -  FuigradW^ 

div(p1t/1/i1  -  Ffygrad/ii)  = 
div(p1L/1/<1  -  F/qgrad/q)  = 

div(p1L?1ej  -  Fsjgradsj)  = 


dP 

=  “ai  +  Ci  ws-wil(wg 

-  W[)  +  rrigiWg 

(142) 

0 

(143) 

Gk\  —  P\0\ 

(144) 

Ci  xrGfci  -  C2pxE\  /k\ 

(145) 

4.2.3.  Interface 

Assuming  direct  condensation  onto  the  surface  of  the  water  jet, 
the  heat  transfer  coefficients  of  direct  condensation  are  given  for 
each  interface  cell  surface,  by  using  the  ideal  condensation  based 
on  the  Clausius-Clapeyron  equation  and  Bernoulli’s  theorem: 


Hci  =  \/2h)g2  /  ((vg(ug  -  u,))(Tsat(Tsat  -  T|i)))1/2  (146) 

mgi  =  Hd(Ts  -  Tn)Ai/hig  (147) 

div(p1L?1/i1  -  Ffygrad/ii )  =  mgl/ig  (148) 

The  interface  velocities  between  the  phases  for  each  interface 
cell  surface  are: 

Wli  =  (rgWg  +  r1W1)/(rg  +  r1)  (149) 


where  rg  is  the  gas  shear  stress  and  X\  is  the  liquid  shear  stress,  C,  is 
the  shear  coefficient  of  interface  between  gas  and  liquid  phases 
and  the  subscript  and  i  is  interface  between  gas  and  liquid  phases. 

They  verified  the  model  by  using  high  pressure  steam  test  data 
(high  pressure  small  size  model),  as  well  as  the  visualized  data 
such  as  temperature  and  pressure  distribution  data.  Good 
agreements  were  found  between  the  test  results  and  the  analysis 
results.  Then  they  used  the  model  to  analysis  the  characteristics  of 
the  large-scale  steam  injector  which  shown  that  the  steam  injector 
driven  jet  pump  could  work  in  the  high  pressure  range  over  7  MPa, 
and  discharged  over  12  MPa  even  at  the  large-scale  rated  flow  rate 
of  220  ton  per  hour. 


The  physical  phenomena  involve  supersonic  flow,  shock 
interactions,  and  turbulent  mixing  of  two  streams  inside  the 
ejector  enclosure  which  are  difficult  to  be  modeled  precisely.  Also, 
the  solving  process  of  the  models  based  on  the  physical 
phenomena  is  complicated  and  time-consuming.  Hence,  in  order 
to  the  evaluate  the  ejector  performance  and  to  instruct  design, 
Huang  et  al.  [43]  derived  empirical  correlations  for  ejector  design 
using  R141b  as  the  working  fluid  after  experimented  on  15  sets  of 
ejectors: 

^-  =  -0.0517^')  +  1.4362  -4.1734  (150) 


A 

=  h0  +  fiirc  +  b2rc2  +  b3rg  +  b4rcrg  +  b5rc2rg  +  b6rg2 

Tint 

+  b7rcrg2  +  h8rc2rg2  (151) 

where  rc  =  P*/Pe ■  rg  =  P\/Pe  and  the  b0  ~  b8  are  coefficients. 

Cizungu  et  al.  [21]  worked  out  following  correlation  based  on 
the  experimental  investigations  of  Dorantes  and  Lallemand  [44]  to 
estimate  the  entrainment  ratio  of  the  ejector: 


co  =  3.32 


1.21\212 

) 


(152) 


where  2.5  <  £  <  4,  \jr  <  6. 

However,  the  empirical  model  cannot  be  extent  to  other  fluid 
and  ejector.  In  order  to  derive  a  relationship  working  for  more 
fluids  and  physical  configurations,  Chou  et  al.  [45]  regressed  a 
relationship  for  the  maximum  flow  ratio  against  the  geometry 
based  on  the  experimental  results  from  literatures  [9,12,46-48]: 


—  =  (^) - 1—nlKpKfKs<p  (153) 

mP  \Pj  (Te/rg)1/2  P 

where  I<p  accounts  for  the  primary  nozzle  performance, 

“’-(rpHri)  >1M> 

wherein  A  is  the  correction  constant  for  the  primary  nozzle 
working  range  and  FPS  is  the  failure  pressure  ratio  of  the  primary 
nozzle.  I<f  is  the  physical  effective  area  ratio  proposed  by  Solokov 
and  Zinger  [49],  I<{  =  fi  x  ff  -  where  ff  =  ^  is  the  ejector 
geometric  design  area  ratio  for  the  ejector,  App  is  the  area  of 
primary  flow  at  choking  location  and  /x  is  a  correction  constant  of 
between  1.35  and  1.5;  I<s  accounts  for  the  “blockage  effect”  of  the 
ejector  configuration,  affecting  the  entrainment  and  mixing 
process  and  cp  is  an  application  specific  correction  factor  which 
can  be  applied  to  an  ejector  operating  with  a  given  refrigerant  and 
system  configuration.  I<s  =  (Ls/A)b~M,  where  Ls/zA  is  the  flow 
aspect  ratio  (Ls  is  the  length  of  the  suction  chamber,  i.e.  convergent 
section  of  mixture)  and  A  is  the  radius  of  A3);  b  is  the  Mach  number 
limitation  constant,  4  for  R113  and  R141b  and  5.2  for  steam,  and  M 
is  the  design  exit  Mach  number  of  the  primary  nozzle. 

More  investigators  have  performed  detailed  experiments  and 
developed  numerous  empirical  correlations  to  predict  the  entrain¬ 
ment  ratio,  gas  hold-up,  mass  transfer  coefficient  and  interfacial 
area  based  on  the  flow  direction:  vertical  up-flow,  vertical  down¬ 
flow  and  horizontal  flow.  Balamurugan  et  al.  [50]  summarized  in 
detail  of  these  correlations  which  will  not  be  covered  here.  As  most 
of  the  models  are  based  on  liquid  as  the  primary  fluid  and  gas  as  the 
entrained  fluid,  Balamurugan  et  al.  [50]  developed  a  semi- 
empirical  model  to  predict  the  liquid  entrainment  rate  (gaseous 
as  primary  fluid)  taking  into  account:  the  compressible  nature  of 
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air  in  the  nozzle,  pressure  drop  for  two-phase  flow  and  the  losses  These  equations  are  valid  over  the  following  ranges:  4500  > 
due  to  change  in  cross  sectional  area.  Ts  >  10°C,3500>Pp  >1001<Pa,  andPexit/Ps  >1.81. 


The  entrainment  ratio  :  co 

=  —7.51  +  9.93 -  2.31  (155) 

The  previously  mentioned  empirical  or  semi-empirical  models 
are  regressed  from  measurements.  They  are  limited  to  the  range 
over  which  it  was  developed,  which  limits  their  use  in  investigat¬ 
ing  the  performance  of  new  ejector  fluids,  designs  or  operating 
conditions.  Hence,  Dessouky  et  al.  [5]  developed  a  semi-empirical 
model  based  on  the  basic  model  similar  as  [10]  and  a  large 
database  extracted  from  several  ejector  manufacturers  and  a 
number  of  experimental  literature  studies  with  attempting  to 
eliminate  the  need  for  iterative  procedures  and  simple  to  use,  also 
give  more  flexibility  in  ejector  design  and  performance  evaluation 
comparing  with  empirical  model. 

The  correlations  for  steam  ejector  performance  evaluation  are: 

co  =  aErbPecPcd  ( e  +  (156) 

(h  +  iPc1) 


oo  =  aErbPecPc 


d(e+  /ln(Pp)) 
(g  +  hln(Pc)) 


(157) 


wherein,  a-i  are  the  coefficients.  The  former  is  used  for  the 
entrainment  ratio  of  choked  flow  or  compression  ratios  above  1.8 
while  the  latter  is  used  for  the  entrainment  ratio  of  un-choked  flow 
with  compression  ratios  below  1.8.  They  also  worked  out 
correlations  used  for  design: 

For  choked  flow: 


Pnp  =  0.13P°33P?73 

(158) 

Am/ Am  =  0.34P,!  09Pp  1 12<w-° 16 

(159) 

Anp/Am  =  1.04Pc-°'83P°86w-°12 

(160) 

For  un-choked  flow: 

Pnp  =  1 .02P-°- 000762  P°" 

(161) 

A1t//4m  =  032P’1Pp1,3«-°36 

(162) 

Anp/AM  =  1.22  p-o-sipo.si^-ao739 

(163) 

They  compared  the  results  of  the  semi-empirical  model  with 
the  experimental  results  derived  from  other  literatures  [6,9,48,51- 
55]  and  good  agreements  were  found. 

Besides,  they  derived  an  empirical  model  for  ejector  perfor¬ 
mance  analysis  [56]  based  on  the  data  and  method  presented  by 
Power  [11]: 


oo  =  0.296 


(Pexit)119 

(Ps)104 


(164) 


PCF  is  the  primary  steam  pressure  correction  factor  and  TCF  is 
the  entrained  vapor  temperature  correction  factor.  They  are 
defined  as: 


PCF  =  3  x  10“7(Pp)2  -  0.0009(Pp)  +  1.6101  (165) 


TCF  =  2  x  10  -  8(TS)2  -  0.0006TS  +  1.0047  (166) 

The  model  equations  are  valid  only  for  an  ejector  operating  with 
steam  as  the  primary  fluid  and  the  entrained  fluid  is  water  vapor. 


6.  Conclusion 

Various  mathematical  models  have  been  proposed  to  assess  the 
performance  of  ejector  for  different  operation  and  design 
conditions.  Generally,  the  mathematical  models  are  based  on 
the  flow  and  mixing  phenomena  inside  the  system.  They  can  be 
divided  into  two  main  categories  in  terms  of  thermodynamic 
model  and  dynamic  model.  Within  the  two  categories,  the  models 
can  be  subdivided  based  on  flowing  phases.  Both  thermodynamic 
model  and  dynamic  model  have  two  types:  single-phase  flow  and 
two-phase  flow.  The  single-phase  flow  thermodynamic  model  can 
be  classified  into  two  sub-categories  according  to  the  mixing 
mechanism:  constant-pressure  mixing  model  and  constant-area 
mixing  model.  The  two-phase  dynamic  model  can  be  further 
subdivided  based  on  the  calculation  methods:  mixture  model  and 
Eulerian  model.  Also,  using  the  measurements  data,  several 
empirical/semi-empirical  models  to  evaluate  the  performance  of 
ejector  and  to  guide  ejector  design  are  reviewed  in  this  paper. 

Thermodynamic  models  are  usually  expressed  in  explicit 
algebraic  equations  and  are  based  on  the  steady-state  one¬ 
dimensional  model,  focused  on  the  pressure  change  caused  by 
supersonic  shock.  In  these  models,  the  detailed  local  interactions 
between  shock  waves  and  boundary  layers,  their  influence  on 
mixing  and  recompression  rate  are  not  taken  into  account.  Lumped 
method  is  applied  in  the  governing  equations  which  lead  to  simple 
model  with  lower  accuracy.  In  order  to  offer  good  accuracy  in 
predicting  the  performance  of  the  ejector  by  the  thermodynamic 
model,  some  isentropic  coefficients  accounting  for  the  friction  loss 
were  applied  in  the  model  and  were  determined  by  experiment.  In 
these  models,  certain  knowledge  on  the  choking,  shock  and  mixing 
should  be  obtained  in  advance  and  assumptions  should  be  made 
for  model  simplification.  The  more  practical  models  known  as 
constant-pressure  mixing  model  and  constant-area  mixing  model 
made  the  model  feasible.  The  constant-pressure  mixing  ejector 
provides  better  performance  than  the  constant-area  mixing 
ejector.  While  the  constant-area  mixing  model  gives  more  accurate 
performance  prediction.  In  order  to  give  more  accurate  prediction, 
the  two-phase  flowing  model  is  developed  considering  the 
condensation  of  the  secondary  flow  or  the  mixture  of  two  streams 
in  different  phase  states.  In  these  models,  quality  is  introduced  in 
the  calculation  of  the  state  parameters.  Though  the  governing 
equations  are  basically  the  same  as  single-phase  flowing  model, 
the  accuracy  is  improved. 

The  dynamic  model,  on  the  contrary,  accounts  the  turbulence 
interaction  between  the  primary  and  secondary  stream,  the  shock 
reflection  and  chocking.  It  is  more  related  to  the  actual  process 
occurred  in  the  ejector  and  the  effect  of  the  geometrical 
parameters  and  operation  parameters  can  be  well  explained. 
The  precision  of  this  model  is  thus  greatly  improved. 

The  flow  in  the  ejector  is  from  supersonic  to  sonic,  and  then  to 
subsonic.  Therefore,  the  choice  of  the  turbulence  model  is 
important.  For  compressible  flowing  model,  it  was  found  RNG  k- 
8  and  k-co-SST  models  were  the  best  suited  to  predict  the  shock 
phase,  strength,  and  the  mean  line  of  pressure  recovery  [31].  For 
incompressible  flow,  the  standard  k- 8  model  and  realizable  k-s 
model  were  widely  used  for  time  saving.  For  two-phase  flowing 
model,  the  mixture  model  can  give  reasonable  results. 

It  can  be  seen  that  the  pressure  inlet  and  pressure  outlet 
boundary  conditions  are  widely  used  in  all  the  models  whenever 
possible.  Other  boundary  conditions  such  as  mass  flow  rate  or 
velocity  are  seldom  used  in  models.  For  thermodynamic  model,  the 
auxiliary  relations  such  as  gas  dynamic  equations,  the  defining  of 
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Mach  number  and  sonic  velocity,  the  state  expression  have  to  be 
employed  to  complete  the  mathematical  model. 

Due  to  the  complex  nature  of  the  partial  differential  equations, 
it  is  necessary  to  solve  the  mathematical  model  using  numerical 
methods.  Finite  difference  method  is  recognized  as  the  most 
accurate  and  most  universal  solution  technique  and  is  widely 
employed  in  ejector  modeling. 

Besides,  model  validation  is  an  important  step  in  model 
development  since  it  offers  the  possibility  of  comparing  simulation 
results  with  actual  system  behavior.  Experiments  are  mostly  used 
to  validate  the  mathematical  model.  Besides,  comparing  with  the 
previous  numerical  results  is  also  a  good  validation  method. 

Though  a  large  amount  of  works  have  been  conducted  on 
modeling  and  analyzing  ejector,  further  efforts  are  still  needed: 

(1)  To  study  the  influence  of  variable  isentropic  coefficients  which 
are  taken  as  constant  in  almost  all  existing  thermodynamic 
models. 

(2)  To  improve  the  accuracy  of  the  model  based  on  turbulence 
model,  since  by  now,  there  still  no  uniform  model  being  used. 

(3)  To  construct  a  simulation  package  of  the  whole  ejector-based 
system  by  combining  the  model  of  the  ejector  and  other 
components  in  the  system. 
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